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Mean-field treatments of Yang-Mills theory face the problem of how to treat the Gauss law 
■ constraint. In this paper we try to face this problem by studying the excited states instead of 
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the ground state. For this purpose we extend the operator approach to the Random Phase Ap- 
proximation (RPA) well-known from nuclear physics and recently also employed in pion physics to 
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• general bosonic theories with a standard kinetic term. We focus especially on conservation laws, 
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^f) . and how they are translated from the full to the approximated theories, demonstrate that the 



operator approach has the same spectrum as the RPA derived from the time-dependent variational 

principle, and give - for Yang-Mills theory - a discussion of the moment of inertia connected to the 
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^ ■ energy contribution of the zero modes to the RPA ground state energy. We also indicate a line of 



thought that might be useful to improve the results of the Random Phase Approximation. 
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I. INTRODUCTION AND MOTIVATION 



Recently, variational calculations based on Gaussian wave functionals have stirred some 
interest when applied to Yang-Mills theory . A problem constantly encoun- 

tered is the fact that Gaussian wave functionals are not gauge invariant, i.e. they are not 
annihilated by the Gauss law operator (except for the case of electrodynamics). Therefore 
one faces the situation, that the mean-field treatment (to which the Gaussian wave func- 
tionals correspond) does not respect a symmetry of the Hamiltonian. This is a situation 
commonly encountered in nuclear physics where one often even wants to break as many 
symmetries as possible in order to store a maximum amount of correlations in a wave func- 
tion of a very simple form. We therefore look to nuclear physics for a possible remedy of 
this problem. 

One possibility is the introduction of a projector which projects the Gaussian wave func- 
tional onto the subspace of gauge invariant functionals |l[ O, S, 0, S, Q . Another possibility 
is the Random Phase Approximation. In this framework one studies not the ground state 
but excited states, and one finds that - if the mean-field ground state has a broken symmetry 
- spurious excitations exist. Under certain conditions these spurious excitations decouple, 
however, and we have the pleasant situation that our excitations are as good as they would 
be if we had started from a gauge invariant mean-field state. 

The computation of the excited states of Yang-Mills theory is especially interesting for a 
variety of reasons: First, one has an alternative point of view from which to shed on light 
onto the confinement problem, compared to the usual approaches that try to compute the 
linear potential between static quarks. Second, one has quite accurate lattice data for some 



This allows 



of the low lying glueballs for both two and three spatial dimensions, e.g. 
to judge the approximations that usually occur in the process of analytical calculations. On 
the other hand, there is quite an amount of information that is very hard to obtain from 
the lattice, like masses of excited glueballs, and also regularities and patterns in the glue- 
ball spectrum which can only be observed but not explained within the lattice framework. 
Third, in recent years the glueball spectrum has also become a matter of interest for theories 
that claim to have a very non-trivial connection to Yang-Mills theory like e.g. supergravity 



theories 



There are basically two approaches to the generalized Random Phase Approximation 
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qRPA). The first one starts from tlie time- dependent variational principle of Dirac 
isl [14 . llS^ . Tfie variational parameters are decomposed into one part that solves the 
static equations of motion and one "small" fluctuating part. In this context, the nature of 
gRPA as a harmonic approximation becomes clear; the general procedure is outlined in app. 
\K\ since the equivalence to the second approach to gRPA has up to now only been shown 
for the fermionic case Il3|. The second approach is based on the formalism of creation and 

n 

annihilation operators, well-known in nuclear physics and has recently been applied 
also in the context of pion physics jl?! . Iisl . Q]. In contrast to these latter investigations, 
we consider here generic bosonic field theories; they are only restricted by the requirements 
that the kinetic energy shall have a standard form, and the remainder of the Hamiltonian 
shall be expressible as a polynomial in the field operators. 

The structure of the paper is as follows: We begin by introducing (in sec. the Hamilto- 
nian formulation of generic bosonic theories which will be the subject of investigation in the 
main part of this paper. In sec. II I II we rewrite the usual canonical formulation in terms of 
creation and annihilation operators thereby bridging the gap between canonical and many- 
body treatment. Especially, we will give an expression of the Hamiltonian for a generic 
bosonic system (with a standard kinetic energy term) in terms of creation/annihilation op- 
erators. This is followed by an alternative form of the Schrodinger equation useful in the 
context of gRPA. We then introduce the two crucial approximations needed to obtain the 
gRPA equations from the Schrodinger equation. The second of these approximations can - 
under certain conditions - be rephrased in terms of the so-called quasi-boson approximation 
which will allow for a simpler formulation. We then discuss the normal mode form of the 
Hamiltonian, and see what kind of difficulties appear if the gRPA equations have zero mode 
solutions. In sec. IIVI we will discuss how conservation laws translate from the full theory 
to the theory approximated by gRPA, when zero mode solutions are implied, and also the 
connection to symmetry breaking in the mean-field treatment is indicated. We will also draw 
attention to the fact that the character of the symmetry under consideration can change 
from a non-Abelian to an Abelian symmetry. In sec. |3 we demonstrate the equivalence 
of the gRPA formulation based on the time-dependent variational principle to the operator 
approach employed in this paper; this equivalence has so far been demonstrated only in the 
case of fermionic systems [l3|. In sec. IVII we discuss further how the normal mode form of 
the Hamiltonian is altered if the gRPA equations have zero mode solutions, and the special 
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role that is played by the moment of inertia. We then restrict the further discussion from 
generic bosonic theories to SU (N) Yang-Mills theory, where we discuss both the moment of 
inertia in general, and also give an explicit computation of the leading terms in a perturba- 
tive expansion. We end the section with a discussion of the possibility of interpreting the 
moment of inertia as the static quark potential. In sec. IVllI we give a critical evaluation of 
the generalized Random Phase Approximation and give an outlook to further applications. 
At the end of the paper a number of appendices is given. In app. El we give a short ac- 
count of the gRPA as derived from the time- dependent variational principle in so far as it 
is needed for the purposes of this paper. In app. El we discuss shortly Yang-Mills theory in 
Weyl gauge, and how it fits into the general framework of bosonic theories as discussed in 
the main body of the text. In app. Oa number of explicit expressions and computations can 
be found, in app. Othe proof for bosonic commutation relations among gRPA excitation 
operators is given, and in app. jElthe proof of a theorem used in sec. Illll can be found. 
The important issue of renormalization will not be addressed in this paper. This topic in the 
context of the gRPA is subject to further investigation. In the approach to the gRPA that 
is based on the time-dependent variational principle a discussion of renormalization issues 



for the (f)^ theory has been given in 



II. HAMILTONIAN FORMALISM 



In this section the Hamiltonian formalism for fairly general bosonic theories is introduced. 
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2l|, we will be very brief here. 



Since there are good references to the subject, e.g. 
The basic variables in the Hamiltonian approach consist of the the field operators (pi and 
the canonical momenta vTj. In order to allow for sufficient generality we will use super- 
indices which can - besides the position coordinate x - also contain spatial indices or internal 
(e.g. color) indices. The Einstein summation convention is adopted, implying sums over all 
discrete and integrals over all continuous variables. These basic variables satisfy the basic 
commutation relations 



yj, — iSij. (1) 

We will work in the Schrodinger representation, i.e. we will not work with abstract states 
but with a field representation of the states (analogous to the position representation in 
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quantum mechanics); the objects to be considered are thus wave functional iIj[4>\, given by 



m = m). (2) 

for the state The states |0) are eigenstates of the field operators (pi. The usage of |0) 
as basis states imphes that we work in the position representation, i.e. in the remainder of 
the paper we will realize multiplicatively and vr as a derivative operator: 

= ct>^m and (0|7r#) = (3) 

I d(pi 

or more pragmatically we will read in all formulas where the field operator (p appears this 
only multiplicative and vr as 6/{i6(j)). 

Since we work in the Schrodinger picture the operators considered are all time- independent. 
The states, however, satisfy the Schrodinger equation 

idtim) = H\m), (4) 

where H is the Hamiltonian of the system. In the main section we will be mostly interested 
in stationary states, i.e. states that can be written as iV'(t)) = e-^^*|V'(0)). These are then 
eigenstates of the Hamiltonian 

H\^P) = E\^P). (5) 



III. FORMULATION OF MANY-BODY LANGUAGE FOR GENERIC BOSONIC 
THEORIES 



A. Creation and Annihilation Operators 



Inthis section we wil 
2 
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_apply and generalize the operator approach used in nuclear physics 
271 l28l . l29l |30| and recently also in pion physics 



We 



draw the connection between the canonical treatment and the many-body language as was 
already briefly indicated in Since for our purposes it is sufficient, and in order to allow 
close comparison to gRPA as obtained from the time-dependent variational principle, cf. 
app. El we will consider Hamiltonians of the form 

1 



H 



+ v[<p]. 



(6) 
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V[(f)] is a functional of the field operator, in the following referred to as 'potential'. 

Since the basis of this approach is the stationary Schrodinger equation, we start therefore 

from the most general time-independent Gaussian state {i,j are super- indices): 

m = Afexp (^-(0 - m\g~' - zS),,(0 - 0), + Z7r,(0 - 0),^ , (7) 

where is a normalization constant. A Gaussian state allows the explicit construction of 
creation and annihilation operators as linear combinations of (pi and vr^ that satisfy the basic 
relations 

aiip[4>] = and [aj,aj] = 6ij, (8) 

where the latter relation fixes the normalization. Using eq. ((7)) as reference state to be 
annihilated by one obtains as explicit expressions for a\ a: 

,t _ TT j l^n-i 



«I = Ui, I \^-Gj,' + 2^S,, J (0 - 0), - z (vr - tt)^. | , (9) 

«^ = Ui, I Qgt,i - 2^S,fc^ (0 - 0), + ^ (vr - 7r)^.| , (10) 

where U is (implicitly) defined via the relation 

UijUjk = Gik (11) 

and could also be called the square root of G. A short excursion on the existence and the 
implicit assumption of reality of U is in order here: Since Gik is a real symmetric matrix 
(matrix is used in the generalized sense s.t. also continuous indices are allowed) one can 
always diagonalize it. Therefore one can also always write down a f/ as given above. However, 
G has to satisfy another condition, namely all of its eigenvalues have to be strictly positive, 
since otherwise one will run into two kinds of problems: If G has a zero eigenvalue G~^ 
does not exist and the matrix element of (tt^) will be infinite. This is certainly undesirable, 
and will be assumed not to be the case in the following. Moreover, if G has a negative 
eigenvalue, is not even normalizable, since in the direction of the negative eigenvalue 
'0[0] will increase exponentially for increasing values of 0. Thus, all eigenvalues of G must 
be strictly positive, therefore U can be chosen to be real. 

Since just given via linear combinations of and vr, one can invert these relations to 

obtain 0, vr in terms of a, and the parameters of the Gaussian wave functional we started 
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with: 

= 0i + f/»j + 4) , (12) 
TTi = TTi + 2i| (^^G;^^ - iJ^ik^Ukja] - (^^G;^^ + zSife)t/fcja^ }. (13) 

One should note that the dependence of the canonical operators on the parameters of the 
wave functional we choose is only seeming, since the creation and annihilation operators 
depend implicitly on these parameters as well. If one inserts eqs. (fTUj) into eqs. (fT^ . 
(|T3|l one obtains an identity = 0, tt = vr. In a practical sense, however, we have transferred 
information that is contained in the wave functional to the operators, since in the following 
the only property of ipi^P] that we will use for practical computations is that aj'?/'[0] = 0. All 
parameter dependence that usually comes about by calculating matrix elements now enters 
the formulas via normal ordering. 

Since we have now a representation of the canonical operators in terms of a, (referred 
to in the following as c/a representation) , all operators permissible in a canonical system 
can be expressed in terms of a,a\ especially the Hamiltonian which for obvious reasons is 
central to the following calculations. 

B. Form of Hamiltonian in the c/a Representation 

We have required the Hamiltonian to have a certain structure [cf. eq. ©]. The Hamil- 
tonian there resolves naturally into a kinetic energy T = ivr^ and a potential term V[(j)] 
which is a functional of only. It is very useful to normal-order these expressions to make 
further progress. In the course of this, one makes the useful observation that one can write 
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the kinetic energy as : 

1 o 



+ 



^7r,2 + ^Tr(G-^)+4Tr(EGS)) | 




'-G^\ + - iU-^m - Um-^) . . ]2a\ a,v 



(14) 



where Tr is a trace over the super-indices. A similar result can be found for the potential part 
(except that it is independent of S and vf). In app. |Elwe will demonstrate that the potential 
can be decomposed into c/a operators s.t. the prefactors can be written as functional 
derivatives of the expectation value of the potential between Gaussian states, and the c/a 
operators always appear in a fixed structure^. In the following we will restrict ourselves 
to the contributions to V with up to four c/a operators, since the terms containing higher 
numbers of c/a operators do not contribute^ to the gRPA matrices that will be introduced 



^ e.g. if in the potential there are terms that contain two creation operators, or two annihilation operators, 
or one creation/ one annihilation operator, they can always be written as factor x (fl\a^j + fliOj + 2a\aj), 
and similarly for all other terms that contain a fixed sum of creation and annihilation operators. For 
the kinetic terms things are a bit different, but we have given the decomposition of the only kinetic term 
allowed in eq. 

^ That this is true can be seen by usage of Wick's theorem and the so-called second gRPA approximation 
a|) = that will be introduced in sec. lIIIDl 
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in eq. 



V[<P] = {V[<P]) 
6 



+ 1 T^^^t*^]) ]Uk,n[^n+<^n 



X [a]^al^a].^ + SaJ^aj^ajg + Sal^aj^aj.^ + aj,aj^aj,j (15) 

By adding the expressions eq. ()14| ITH|l together, one observes that the Hamiltonian has a 
very simple schematic structure: 

H = {H) 

+ I -^7^ + - - UJ:U-^) . .]2a] a,, 

+ Q I ^^"^ ) ^kinUk2j2Uk3j3 

+ fi (77^ A?^ ^kinUk^jiUkajsUk^j^ (16) 
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One should note that in all save the last term, the derivatives are taken of (H) whereas in the 
last term the derivative is taken of {V). The importance of this will become clear in sec. lVBl 
Up to now, we have taken an arbitrary Gaussian as a reference state: we have defined our 
creation/annihilation operators relative to that state - nothing else. Eq. is exactly the 
same as H given in eq. (jHI). However, it is obvious from eq. that the Hamiltonian greatly 
simplifies if we choose a specific reference state: a state that is a stationary point of the 
energy functional (H) under variation of the parameters 0, tt, G, S - in other words, a state 
that satisfies the Rayleigh-Ritz variational principle^. 

C. Alternative Form of Schrodinger equation 

In the following we want to generalize the RPA known for many-body physics to the 
quantum field theory described by the Hamiltonian of eq. (jH)). For this purpose we assume 
that the exact (excited) state {u) can be created from the exact vacuum |0) by an operator 
Ql, i.e. 

W)=Ql\o). (17) 

If El, denotes the corresponding eigenvalue of H we have 

HQI\0)=EM\0). (18) 

If we denote the vacuum energy by Eq, we can write this equivalently with a commutator: 

[H,Ql]\0) = iE,-Eo)Ql\0). (19) 

We may multiply both sides of the equation with an arbitrary operator 6Q, and obtain an 
expectation value by multiplying from the left with (0|: 

{0\6Q[H,Ql]\0) = {E,-Eo){0\6QQl\0). (20) 

By subtracting zero on both sides we obtain an equation that only contains expectation 
values of commutators: 

{0\[6Q,[H,Qlm = {E,-Eo){0\[6Q,Ql]\0). (21) 

^ One should note that we require the same in the time-dependent approach, when we decompose the 
time-dependent parameters 7r(i), G{t), T,{t) into a static part and small fluctuations, cf. app. ^ The 
equations for the static part are identical to the Rayleigh-Ritz equations. 
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Using this equation as a starting point to derive the (generahzed) RPA is known in nuclear 

nn 

physics as the equations of motion method |16|, |32| . 



D. First and Second gRPA Approximations 

The generahzed Random Phase Approximation consists now of approximating eq. 1)2111 
in order to obtain a solvable set of equations. Two obvious candidates for approximations 
suggest themselves: first, the excitation operator Ql, second, the vacuum state. In the re- 
mainder of the paper, we will call the approximation concerning the first topic the first gRPA 
approximation, the approximation concerning the vacuum state the second gRPA approxi- 
mation (even though this may be doubtful linguistically). The second gRPA approximation 
is actually easier to state, therefore we start with it: approximate in all expressions involving 
vacuum expectation values of commutators the true vacuum state by the reference state of the 
creation/annihilation operators, which in this context we will usually refer to as mean- field 
vacuum and which will be denoted by | ) . Later on, we will be even more restrictive and 
require the reference state to really be a stationary state of the Rayleigh-Ritz principle as 
was indicated before, but we will state explicitly from when on this additional restriction 
will be necessary. As already noted, the first gRPA approximation deals with the class of 
allowed operators. In nuclear physics it is quite reasonable to assume that the lowest excited 
state above a Hartree-Fock ground state consists of a particle-hole excitation. This results 
in the so-called Tamm-Dancoff approximation. In the generalized RPA one assumes that 
one has a correlated ground state, s.t. not only the creation of a hole and a particle, but 
also their destruction is a possible excitation. In Yang-Mills theory - the theory we will be 
ultimately interested in - it is by far not so clear what structure the lowest excitation will 
have; we have thus taken the two following guiding principles (an argument similar to our 
second principle can be found in jl?! . loj ]) 

1. one of the main differences between fermionic and bosonic systems is that in the latter 
there exist single-particle condensates; thus, one has at least to extend the ansatz for 
the excitation operator by linear terms allowing for fluctuations. Furthermore, since 
we can compare to the gRPA as derived from the time-dependent variational principle, 
we will see that the ansatz for Ql to be proposed leads to equations of motion that 
are identical to those derived from the time-dependent variational principle, thereby 
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verifying the ansatz to be the correct one. 

2. The principal goal of the Random Phase Approximation is the restoration of symme- 
tries that are violated at the mean-field level. We give here a short outline, drawing 
from concepts that will be introduced further below, in order to motivate from this 
property of symmetry restoration the form to be allowed for the excitation operators. 
The first point is that the second gRPA approximation can be replaced (under certain 
conditions) by the so-called quasi-boson approximation (QBA). There, it will turn out 
that gRPA equations can be written as 

[Hb,Q^Bu] = {E.-Eo)QL (22) 

where Ob indicates that the QBA has been used. One the other hand, for a sym- 
metry generator (in our case the Gauss law operator) P which is a one-body oper- 
ator (i.e. an operator that can be written as a linear combination of the operators 
a, a^, aa, a^a\ a^a), one can derive that 

[H,T] = 0-^[Hb,Tb]=0. (23) 

If the mean- field vacuum is annihilated by P this is a trivial statement, since in this 
case Tb = 0. However, if the symmetry is violated on the mean-field level, one obtains 
a non-trivial result, namely that Pg is a solution of the gRPA equations with zero 
excitation energy provided that the class of excitation operators contains Tb- Since 
all gRPA excitations are orthogonal this would be a very desirable state of affairs: 
effectively the spurious excitations caused by the symmetry violation on the mean- 
field level would not affect the physical excitations we are interested in. 

To put it all in a nut-shell: one looks at the Gauss law operator, determines its structure (cf. 
app. IG 2|) . and keeps the class of excitation operators so large that the Gauss law operator 
belongs to this class. All the concepts and claims will become clear and will be proved in 
what follows. After this rather long motivation, we will consider excitation operators of the 
following form: 

Ql = l {XL<^14 - Y:^,a^a^ + {Zi:^al - Z'^a^) , (24) 

where X, Y, Z, Z are called amplitudes. From the form of the Gauss law operator one would 
have expected a few more terms; the reasons for leaving them out will be explained in 
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sec. IIIIFI With the ansatz (H^i we derive from eq. (PT|) a closed set of equations for the 
amphtudes X'^, V^, Z", Z". To this end, the (up to now arbitrary) operators 6Q will be 
appropriately chosen, together with the second gRPA approximation. We choose 6Q in such 
a way that the different amplitudes are extracted individually on the RHS of eq. (j2H) . It will 
be useful to consider 6Q G {oj U, and compute the commutators [SQ,Ql], 

approximating the real vacuum by the mean-field vacuum according to the second gRPA 
approximation. We obtain: 

([ataj,Qt]) = (Ka„Qt]) = (25) 

where we have used the notation (| . . . |) = (...) and introduced the abbreviation 
defined as Xjjj} = Xij + X^j, and correspondingly for Y{ij^- With the bQ^ defined above 
we now have to compute the LHS of eq. 1)211) . For this purpose we introduce a number of 
matrices 

^nimi = {[ala],[H,al^al]]), Dnjm = {[ala],[H,aJi]), 
Bnjmi = {[aia],[H,amai\]), E^m = {[ai,[H,alJ\), (26) 
Cnjm = {[ala],[H,al^]]), F„„ = ([4, [iJ, aj]), 
and study their properties under interchange of labels. This will allow to reduce the num- 
ber of independent entries of the LHS of eq. ()2H). The basic tools for this study are the 
Jacobi identity and the second gRPA approximation, i.e. a\) = 0. We also use frequently 
that [a, a] = [a^,a^ = 0. It turns out that the matrices A and E are symmetric in all 
indices, whereas F is a hermitean matrix, F*^ = Fmn- In order to establish that B is 
also hermitean in the sense B'^j^^ = Bminj (whereas it is symmetric under interchange 
n ^ j,m ^ i) we need the second gRPA approximation since only upon usage of this 
approximation {[H, [aj„a|, a„aj]]) will be zero generally. 

Apart from one matrix, all other matrices that arise from inserting the different 6Qs into 
the LHS of eq. ()2ip are trivially related to the matrices A, . . . ,F introduced above; this one 
non-trivial matrix is 

{[al[H,ama,]]). (27) 
It is related to the matrix D, cf. eq. ()26p . via 

{[al[H,a^ai]]) = 6inAm + Snm^i + D*^.^ , (28) 
13 



with Aj^m = {[H, ai^m])- In sec. lIII El the second gRPA approximation will be replaced by the 
so-called quasi-hoson approximation. In that formulation, it will be obvious that we have to 
neglect the terms A. From the explicit expression for the expectation values, given further 
below in eq. (j38|) . we see that, at the stationary mean-field point^, A is indeed zero. We note 
two points: 

• The gRPA matrix [essentially the LHS of eq. ^2,\\ \ is hermitian iff A = 0. 

• The second gRPA approximation and the quasi-boson approximation to be introduced 
below are compatible iff A = 0; such a constraint does not appear in nuclear physics as 
is easily comprehensible in two different ways: First, the matrix under consideration 
(as well as the matrices C, D) are non-zero only if the Hamiltonian contains terms with 
in total three creation/annihilation operators, i.e. terms that violate particle-number 
conservation. Such terms are not allowed in the usual nuclear physics framework, and 
thus there is no possibility for the above consistency condition to arise^. Second, if 
we replace (as we will do in sec. IIIIE|) the operators containing two ordinary boson 
creation operators by a new boson operator, in the case of bosonic theories we have 
still to keep the original boson, in contrast to the RPA treatment of fermion systems, 
where it is not necessary to retain the original fermions once one has the formulation 
in terms of bosonic operators at hand; 



This whole discussion allows us now to write down a first form of the gRPA equations: 
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with VLjj = Ey — Eq. The factors | associated with some of the matrices C, D shouldn't 
lead one to conclude that the gRPA matrix on the LHS is not hermitian. In fact, as we 
will see later on, we can write the gRPA equations in a very compact form resulting from a 



i.e., that point in parameter space where the energy expectation value is stationary with respect to 
variations of ^, tt, G, S; for A to be zero it is necessary that -^{H) — 0, and j4-(i7) = 0. 
^ In the bosonic theories under consideration here, these terms with three creation/annihilation operators 
appear in the Hamiltonian usually due to a condensate. 
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hermitian Hamiltonian (at the stationary point of the mean-field equations -~{H) = 0, and 

s 



E. The Quasi-Boson Approximation 

Sometimes, it is useful to have the second gRPA approximation at hand without having 
to take expectation values. This can be done with the help of the quasi-hoson approximation. 
The name has its roots in nuclear physics Q, 2^, where two fermions are combined into 
one boson. This works only approximately^, so that the result was christened a quasi-hoson. 
Interestingly enough, we can do the same in a bosonic system: we replace a two boson 
operator (like aa or a^a"!") by a new operator that again has bosonic commutation relations. 
We construct the new boson-pair operators B, B"'" s.t. the commutation relations are identical 
to the mean-field expectation values of commutators containing still the pair of original boson 
operators: 

[amai.aj] =0 [B™,aj] = 

[amai^a]] = 6ijam + Smjai ([amai,aj]) = ^ [Bmj,aj] = 

[a^alaj] = -6ija\^- 6.mja\ {[a^^a\,aj\) = ^ [BL^^j] = 0. 

The only non-vanishing commutator involving B, B"!" originates from [0^0^, aj^aj]: 



^rrw; •-'nj 

We therefore replace 



^(a-mO-i) ^ Bmi ^^^^ 



The operators B, B^^ have the following commutation relations 



[aj, Bmi] — [aj, Bj„j] — [B^j, B„j] — (31) 
and [B^i,By = ^(5m5„j + (5i/^„). (32) 



^ i.e. the ordinary bosonic commutation relations are only fulfilled if we take the mean-field expectation 
value of the commutator of the boson pairs, and thus in some sense employ the second gRPA approxima- 
tion, by using the mean-field instead of the exact vacuum state. 
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With these new operators, the excitation operators 
become: 

At this point a short comment on how one constructs quasi-boson approximations is in order, 
cf. also eq. below. Usually we will have to consider only one- and two-body operators, 
but the procedure should work for higher-body operators as well: 

(i) first, one writes down the operator O under consideration as a polynomial of a, 
operators. Then one tries to extract the coefficients of the different powers in analogy 
to eqs. (j2Sl), via (a) taking multiple commutators of O with up to two a, a) oper- 
ators, and then (b) taking vacuum expectation values of these multiple commutators 
employing the second gRPA approximation. 

(ii) The quasi-boson approximation Ob of this operator O is then given as a polynomial of 
B, ^\a,a^ operators. The coefficients of the different powers are extracted by taking 
multiple commutators of Ob with a, a) , B, B"'" without taking expectation values. One 
then requires that the coefficients oi Ob determined in this way are identical to those 
of O determined in (i) if one replaces in the multiple commutators of (i) O hj Ob and 
terms of the structure aa by ^/2B and a^a^ by -\/2B"'' as indicated in eq. 

This is precisely the procedure that led to eq. Using the procedure outhned, one can in 
fact express the QBA Ob of a general one-body operator O in terms of mean- field vacuum 
expectation values of commutators: 

Ob = {0) - l={[ala],0])B,, + ^{[a,a,,0])Bl^ - {[alO])a, + {[^^^ (35) 

Incidentally, one can at this point easily verify -by writing out the commutators - the claim 
made above: the QBA of a symmetry generator (which is generally hermitian) of one-body 
type that annihilates the mean-field vacuum vanishes. 

Now we will follow the above procedure to obtain the Hamiltonian in quasi-boson approx- 
imation. The coefficients mentioned are nothing but the matrices A,...,F. We therefore 
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have the requirements 



A. 



D. 




njm 



nm 



'nm 



(36) 



where Hb denotes the Hamiltonian in quasi-boson approximation. At this point it seems in 
order to discuss a question that turned up in connection with eq. ()27|) : There we saw that 
one of the double commutators was not directly connected to one of the matrices A, ...,F but 
that some extra terms appeared, the 'A-terms'. If we consider the same double commutator 
in the quasi-boson approximation, we obtain 



=0 

Thus the quasi-boson approximation implies the vanishing of A. An explicit expression can 
be given for A as well (cf. app. IC 1|) 



which vanishes if we choose the parameters of the reference state s.t. the energy is stationary 
under their variation (at least w.r.t. 0, vf). We obtain a consistency condition between the 
second gRPA approximation and the quasi-boson approximation: the latter is equivalent 
to the former only at the stationary point - at least w.r.t the condensates 0, tt - of (if). 
Thus, in the following we will always assume that the state (which we often call the mean- 
field vacuum) relative to which our creation/annihilation operators are defined [cf. eq. (j?))] 
has its parameters chosen s.t. it satisfies the Rayleigh-Ritz principle - that this restriction 
would become necessary was already indicated when the second gRPA approximation was 
introduced. With this qualification, we can give the (now hermitian) Hamiltonian in quasi- 
boson approximation: 

ff p 1 (a d d _i_ 4* I _i_ 1 R R 




(37) 




(38) 




(39) 



2 {Enin2'^niO'n2 + -^nin2'^rii '^^2) ~^ -^»t2ril '^ni 
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with Emf = (H). We have achieved writing the approximated Hamiltonian as a quadratic 
form, which can always be diagonahzed. Using this Hamiltonian Hb one can see that the 
gRPA equations eq. (|29|) can be written in the transparent form^ 

[Hb,QU = ^^QL- (40) 

This form of the gRPA equations demonstrates that we have performed the approximation 
to the dynamics (Hamiltonian) and to the excitation operators consistently, since the form 
of the equations is identical to the form of the Schrodinger equation eq. ()19|) with the exact 
Hamiltonian and the exact excitation operators. 



F. Restriction on Excitation Operators, part II 

Having introduced all the concepts of gRPA, it is time to look again at some terms that 
appear neither in Qj^ nor in Q^^. 

1. Ql does not contain an a^a term. The superficial reason for this is that the correspond- 
ing amplitude cannot be extracted in a way similar to the other amplitudes X, Y, Z, Z, 
since due to the second gRPA approximation 

{\[6Q,a^a]\) = {\6Qa^ a\) - {\a^ a6Q\) = 0. (41) 

=0 =0 

This point will become clearer in the context of the QBA, cf. item IHl below. 

2. Ql does not contain an aa^ term. This is due to the fact that aa^ is distinguished from 
a^a only by a constant. However, this constant would lead to a non-vanishing expecta- 
tion value of Ql between mean-field vacua, in view of the second gRPA approximation 
a situation certainly not desirable for an excitation operator (all other terms contained 
in Ql vanish between mean-field vacua !). One could also argue on a formal level that 
the way the generalized RPA is derived here does not allow to determine any constant 
parts of Ql, so if we cannot determine a constant, we shouldn't put it into our ansatz 
in the first place. 

''' One inserts eq. H39() and eq. (|34|l into eq. (|40|l and compares the coefBcients of the different operators 
a, a'f, B, B'l^. This then reproduces eq. (f^ . 
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3. Here we want to have a look at the fact that o)a does not appear in our excitation 
operator from the perspective of the QBA. The excitation operator is hnear in all the 
different boson creation/annihilation operators that we construct as one can see in 
eq. (j34p . Thus, the question arises whether we can construct a boson operator from 
a^a similar to how we construct one from a^a) . The problem is immediately apparent: 
if we consider Cmi = ctl^CLi, then its adjoint has the same structure as Cmi, since 
^iii — (^I'^m = Cim] if we compute the commutator [Cij,Cli] = 5jka\ai — 6iia\aj we 
see that its mean-field expectation value is zero. Therefore we cannot construct boson 
operators with the correct commutation relations from a^a in the same way as we did 
for o)a\ and thus they don't appear in the excitation operators Q\. In nuclear physics 
this is well-known. There one can see that the quasi-boson operators correspond to a 
creation of a particle-hole pair (or its annihilation), whereas the operator that creates 
and annihilates a particle is translated into an operator that creates and annihilates 
a quasi-boson {2^ . 

G. Normal Mode Form 

In app. Owe consider the commutation relations of the normal mode operators. There 
we show that if all excitation energies Vty are distinct and non-zero, the normal modes have 
the usual bosonic commutation relations: 

[Qbm,QL] = '^/- (42) 
Using this fact, we can conclude from eq. (pUj) 

Hb = Erpa + J2 ^-QLQbu. (43) 

V 

The sum over v extends over all positive semi-definite The constant Erpa can be 
determined as usual, cf. namely by requiring 

(Hb) = Emf. (44) 

Using 
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with yf„,} = + this gives 

erpa = emf - ^ E i^^wI' + E ■ (46) 

fc" \ mi i / 

A similar expression can be obtained in nuclear physics j3| but without the appearance 
of l^j^p. The changes brought about by excitations with zero excitation energy will be 
considered in sec. IVII 



IV. CONSERVATION LAWS 



In this section we want to demonstrate that - under certain conditions - conservation 
laws from the full theory translate to conservation laws in the gRPA treatment. The proof 

n 

follows the lines of the fermionic case as given in |22;]. For the conservation laws to hold in 
the approximated theory, it is mandatory that the reference state minimizes the expectation 
value of the Hamiltonian, since then the terms of the Hamiltonian linear in a, a\ aa, a^a^ 
vanish, cf. eq. One should also keep in mind that the terms of H linear in a,a^ have 
to vanish anyway, since otherwise a treatment using the QBA is not valid as was discussed 
before. 

The fact that conservation laws translate into the approximated theory has to be taken 
with a grain of salt, however, since it turns out (and will be discussed in sec. IIVB|) that 
the character of symmetries may change. Non-Abelian symmetries are usually reduced to 
Abelian symmetries which is not surprising since the gRPA is basically a small-fluctuation 
approximation which does not probe the group manifold. 

A. General Observation 

With these qualifications, let us now compute the commutator of Hb, cf. eq. fjHnjl . with 
a general one-body operator Tb 

Tb = To+ {Tl^A. + r°>„J + V2 {rZnAn, + Cn.B„,„,) , (47) 

where the coefficients have been chosen, s.t. the full operator before QBA reads 

r = rn + r^°a"^ +r°^a +1^° a'^ +r°2 a a +r^^ a (48) 

^ J- ni^Tii ' ^ ni"ni ' ^ ni7i2"'ni"'n2 ' ^ 711^2 ""i ""2 ' ^ nin2"ni""2- \^°J 
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im2 



Then the commutator gives 

[Hb, Tb] = ——j=Bniji {Cnijimi^mi + ^nijimi^mi ~^ ^nijimim2^ mim2 ~^ ^nijimim2^ mim2) 

I ^ pt ( r^* F*^^ -I- Pi* F^*^ -I- A* F"^ -I- R* F^'^ 

' ^^y^^ niji\^nijimi mi ' nij\m\ mi ' nijimim2 mim2 ' nijimim2 mim2 ) 

—n (W F^O -I- P fO^ -I- F^O r" -I- f02 n* 

"jii l^-C/n.imiJ- mi ^riimi'- mi mim2 ^mi"i2ni mim2^mim2ni) 

-A-n"^ ( P* F"^ -I- F* F^*^ -I- F*^^ (^* 4- V^^ D \ (AQ\ 

'"'ni \-^nimi^ mi ' "'^nimi mi ' mim2^mim2ni mim2 'ni"i2".l ) ■ \^^) 

Let us now concentrate on the first hue to make the principle clear; we use the definition 
of A, B, C, D as mean-field expectation values of double commutators involving the full 
Hamiltonian, cf. eq. (IS^ . Then it is clear that - using also eq. (j^Hj) - one can rewrite it as 

Bni,i [H,al^]])rZ + {[ala]^, [if,ami]])C 

+ ([ani4l' (^lni4n2Wmim2 + ([4i4l' t^' ]] )r^^,, 

= B„,,, (([4^4^, [if, r]]) - T'l^^ {[ala]^, [H, al^a^,]]) ) . (50) 

^ . ' 

(*) 

The last term (*) is proportional to i.e. that term in the Hamiltonian that multiplies 
two annihilation operators, and is zero in the case we are considering here, i.e. our reference 
state minimizes the expectation value of the Hamiltonian. The procedure can be repeated 
for the other three terms in [Hb,^b]- Eventually, one finds 

[Hb^Ts] = -^{[ala]^, [^,r]])B„,,, + -^(Ka,„ [H,T]])BI^^ (51) 
— {[al__^, [H, r]])a„^ + {[ttni, [H, r]])a]^^ + terms that vanish at the stationary point. 

Two points should be noted: 

1. In eq. ()H5j) we have given a formula of the QBA of a general one-body operator; if 
[H, r] is a one-body operator, then by comparing eq. ^T\i to eq. (jH^ we conclude that 
[Hb,^b] is the QBA of this operator (apart from a possible mean-field expectation 
value). 

2. The stationarity condition [H'^'^ = H^^ = 0) has been used several times; if we are 
not at the stationary point, some of the terms neglected above do not vanish, i.e. 
[H, r] = does not imply [Hb, F^] = away from the stationary point. 
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We want to assume that we have chosen the reference state s.t. indeed H"^^ = H^^ = 0. 
Then eq. f|FI| simphfies to 

[HB,rB] = 0. (52) 

Comparing this with the gRPA equations as given in eq . (j4m) we see that - since Tb is exactly 
of the form of Q^^^ - that is a solution of the gRPA equations with excitation energy 
zero, i.e. it is a zero mode. 

B. Change of Symmetry Character 

It seems that by this construction we have precisely obtained what we wanted: the 
spurious excitation made possible by the deformed mean-field state is a solution by itself, 
and does not influence the other - physical - solutions inappropriately. However, this result 
has to be taken with a bit of caution: as has been indicated before, the character of the 
symmetry may change. This can be seen as follows: Compute the commutator of the QBA 
Qb,Pb of two arbitrary one-body operators Q, P. Using eq. it is very simple to obtain 

[Qb,Pb] = {[Q,P]). (53) 
In the case of Yang-Mills theory, we obtain for the commutator of the Gauss law operators 

[r^,r^] = ([r,r^]) = ^r''^(n. (54) 

At the stationary point, however, we have for the bosonic theories under consideration 
7f = S = 0, since (vrvr) is quadratic in both vf and S, cf. eq. ()A7|) . Thus, in the case of 
Yang-Mills theory, we obtain 

[r^,r^]=o, (55) 

using the expression given in eq. ()C16|) for (F'*). In other words, the QBA has reduced the 
non-Abelian SU{N) symmetry to an Abelian f/(l)^^~^ symmetry. This is a phenomenon 
also well known from perturbation theory. 

C. Possibilities for Improvement 

One main drawback of the gRPA is the fact that non-Abelian symmetries may be reduced 
to Abelian symmetries. There have been a couple of investigations in nuclear physics, e.g. 
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ly,|2^, and other quantum field theories (where the fields were in the fundamental repre- 







18l igj l whether one can regard the gRPA as a certain order in a systematic 



sentation) 
expansion^, symbolically 

complete result = mean field + gRPA + higher orders. (56) 

A first hint on how one may construct such higher orders for one-body operators like the 
Gauss law operator can be obtained from the following observation: if we calculate (as we 
have done above) the commutator of two quasi-boson approximated Gauss law operators, 
we obtain the mean-field value of the full commutator; apparently, in the expansion alluded 
to above, taking a commutator reduces the order in the expansion by one (for a similar 



observation in nuclear physics, cf. e.g. 



Ifl |3|). If we construct a next order of say 



the Gauss law operator, taking the appropriate commutator should give the quasi-boson 
approximation of the commutator, the latter being again a one-body operator. In formulas, 
if we write^ 

^_r^ = (r) + Tb + T2B +higher orders (57) 

full operator mean-field conventional gRPA next term 

we have 

n,r%] = {[r-,r>']) (58) 

and T2B would be correct if it would reproduce 

[rB,r^B] + [r2B,r^] = r,r%. (59) 

Such a T2B can indeed be constructed. If one requires that T2B reproduces the same com- 
mutators as Hb in eq. (jH(3|) [with H obviously replaced by T on the LHS of eq. one 
obtains T2B as 

r2B = 2ri]BLB,, + Tf^,,alal + r°5,a,,a,, + rga^a, (60) 
where the coefficients F^j* are defined in eq. (jIH|) . 

For the calculation of the commutator [F^, F2^] we use a technique very similar to the one 



® Within the path integral approach one can show that the mean-field approximation and the RPA corre- 
spond to the first and second, respectively, order in the loop expansion |33|. 

^ One should note that this is a slight deviation from the conventions in the rest of the paper, e.g. in eq. (|35|l 
the mean- field order (O) is assigned to Ob- 
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used in sec. II V Al by rewriting the coefficients of T^^ as double commutators (note that the 
coefficients of T^^ are called 'Tj^ etc; note also that - as in the remainder of the paper - we 
use T{ij} = Tij + Tji) 

Using these results, one can rewrite 



= [r^,r^^] (62a) 
-a,X[at^,[^^n])+at^(K,,[^^n]) (62b) 
-i=B,,,,([a;a;, [r,r^]]) + ^Bt^^.^(K,a,„ [^^^^]]) (62c) 

Since [r'^,r*] just again gives a one-body operator, we can compare eqs. ()62b|) . ()62cp to 
eq. dSni), and see that it is indeed the QBA of^° [T", T''], i.e. [T", r%. Thus we have fuffiUed 
our goal: we have constructed an extension of r% which improves the commutation relations 
as was desired. 

At this point one should keep in mind the following: We have constructed by requiring 
that ([ci|aj, [r,afc]]) = [-\/2BJj., [r2_B,afc]]. This gave zero as the coefficient of a'B. Had we 
chosen instead to require ([a^, [F, a|at]]) = [a^, [T2B, V^^lj]] we would have obtained a non- 
zero coefficient for a^B. This incompatibility of the two construction principles goes back 
basically to the problem that the commutation relations for a, B are constructed to work with 
one-body operators (as was noted in the context of Hb where both coefficients are required 
to be equal. This can be achieved by going to the stationary point since the offending terms 
all vanish there; here we cannot follow the same route since the state is already specified). 
The ambiguity is here in some sense resolved, since the latter construction does not lead to 
the aspired commutation relations^^. 



Note again that in this section the expectation value ([r°,r^]) is not incorporated into [r°,r'']B. 
There can be circumstances where even in the latter construction the correct commutation relations can 
be produced. We have investigated the case where the symmetry generators are the Gauss law operators, 
and found two instances where the correct commutation relations were produced even in the presence of 
the aB^, a^B-tcrms. These examples were: G^^ cx 1 - this is not a very likely mean-field vacuum since 
there is no correlation between different points in space - and the case where Df7~^ — 0. For the definition 



24 



V. P-Q FORMALISM 



In sec. IIII Gl we have learned that the usual gRPA procedure encounters problems when 
the gRPA equations have zero mode solutions; in sec. IIVI we have learned that in the 
case we are interested in - namely a mean-field vacuum that breaks certain symmetries - 
zero modes appear necessarily. Therefore, it seems that we may have a problem here. This 
problem, however, was already encountered in nuclear physics, and solved in 3]- We present 
the problem and its solution for the simplest case available, namely the simple harmonic 
oscillator. We then give an alternative formulation of the gRPA in QBA using - instead 
of creation/annihilation operators - operators with canonical commutation relations. This 
formulation will furthermore allow a simple comparison between the gRPA equations derived 
from the time-dependent variational principle and the operator approach. 



A. Resolution of Zero-Mode Problem 



In the nuclear physics literature it is well-known 



that the appearance of zero-modes. 



despite being welcome in our case, leads to all sorts of problems, especially that the solutions 
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of the gRPA equations no longer form a complete set^^. Marshallek et al. 
pointed out that the problems one has are actually an artifact of the creation/annihilation 
operator formalism. We want to give a very primitive illustration of their point: One can 
write down a quantum mechanical oscillator in two equivalent forms (we have set the mass 



of C/, cf. eq. But this requirement on U^^ implies that G^^ has zero modes, and is thus excluded 

in the gRPA framework, cf. sec. IIII Al Thus, for practical purposes, at least in Yang-Mills theory, the 
required commutation relations resolve all ambiguities in the construction of T2b- 

In nuclear physics this is a problem since one there cannot prove that the Hamiltonian can be written as 
a sum of Q\jQu- Completeness is lost in the following way: We see - in eq. l|Dip - that from every solution 
of the gRPA equations we can form another solution by considering the adjoint operator, therefore the 
solutions to the gRPA equations always come in pairs. In the case of zero modes, this is not necessarily 
true any more, since now the 'excitation operators' that we find may be hermitian; thus, by considering 
the adjoint, we do not get a new solution. In fact this situation is the usual case for zero modes, as we 
have seen in sec. IIVI the excitation operators for the zero modes are just the symmetry generators, which 
are usually hermitian. In our discussion in app.^ we hit upon the problem from a different direction; we 
see that only if there are no zero modes in the spectrum we can show that the excitation operators can 
be considered as independent oscillators. 
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equal to one for simplicity): 

^ = y + T^' (63) 

= ^(ata + i), (64) 

where uj denotes the frequency of the harmonic oscillator. If we now send the frequency to 
zero, something seemingly strange happens: 

uj{a^a + i) 0. (66) 

If one now looks at the transformation from canonical to c/a operators, one begins to see 
more clearly: 



and if we take this small-u; behavior into account, there is actually no paradox , since 

—a' a = —[—=)p = — . (69) 



2 2 
This is actually the solution of the problems involving zero modes: one puts the system into 
an external field, s.t. there are no zero modes any more, since the external field^^ breaks 



One should note that, in this hmit ^ 0, no new information can be obtained by considering the adjoint 
of at. 

^'^ The external field can in fact be engineered just to do that, and this works as follows: We solve the 
gRPA equations and determine all modes available. The operator belonging to the zero mode will be 
known beforehand - it is just the quasi-boson approximation of the (known) symmetry generator, which 
is usually a hermitian operator, and so will be the operator associated with the zero mode. One then 
constructs a canonically conjugate variable, s.t. [q,p\ = i and q commutes with all other normal mode 
creation/annihilation operators, and adds a term w^q^ to the Hamiltonian; this will lift the zero mode 
to a finite frequency. By this construction, only the zero mode will be lifted and nothing strange will 
happen to the other modes. Two points should be mentioned: First, the construction of a canonically 
conjugate coordinate to a such a generator is only possible in the gRPA context, since we have seen in 
sec. II VI that the symmetry we have in mind is usually reduced to an Abelian symmetry; if we would still 
have the non- Abelian symmetry this would not be possible (the problems with the construction of 'angle 
operators' in nuclear physics are well-known, for a recent discussion cf. e.g. [sfij). As a second point, one 
should mention that in a quantum field theory there might be problems with breaking gauge symmetries 
in intermediate steps. Since we have indicated here only a formal construction, it is not easy to see where 
these kinds of problems might surface. 
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all symmetries. One diagonalizes the Hamiltonian into a set of harmonic oscillators using 
canonical coordinates, and at the end sends the external field to zero. Then the would-be 
zero modes indeed become zero modes with the correct operator representation, namely p^, 
and the non-zero modes can be rewritten using creation and annihilation operators. As a 
last comment on this topic, one should keep in mind the following: since we already have 
some prejudice about the zero mode operators (as the symmetry generators), we cannot just 
choose their mass to be one as we did in the example. We will rather have to allow for a 
mass tensor, and this will be dealt with in section IVII 



B. Equivalence of Operator and TDVP Approach to gRPA 

We have seen that at least three formulations (using the tdvp, the two gRPA approxima- 
tions, and the quasi-boson approximation) of the generalized Random Phase Approximation 
of bosonic systems are available. In this paragraph we want to add one more, since it sim- 
plifies the treatment of those Hamiltonians we set out to consider, cf. eq. Since we start 
from the gRPA formulation using quasi-boson operators, the whole treatment is only valid 
for the mean-field vacuum representing a solution of the Rayleigh-Ritz variational principle 
(at least w.r.t. variations of (f) and vf). In this section we assume that the mean-field vacuum 
is a solution of the Rayleigh-Ritz variational principle w.r.t. all parameters 0,7f, We 
consider linear combinations of a, a) and B, s.t. 



Qmi — '~f^{^mi ~^ ^mi) Pmi — '^{^mi " ^mi) ■ 

They fulfill the usual canonical commutation relations: 



(70) 



2 

and all other commutators vanishing. If we now rewrite the Hamiltonian Hb in terms of 
these new operators, we obtain a complicated expression where the real and imaginary parts 
of the matrices A, ...,F are separated. We can now put to use that we have computed the 
matrices A, ...,F in app. K^ II ; we have found that 

1. A,C,D are always real 

2. C = -D 
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3. B, E,F may be imaginary; however, as this imaginary part is due to S, and in the 
theories under consideration^^ at the stationary point S = 0, the imaginary parts of 
B, E, F vanish. 



4. As a matter of fact, E vanishes altogether at the stationary point. 

5. A last simplification is that one may decompose B = B^^+B"^"^ and, using this notation, 
we have A = — i?^^. 

Using this information, the Hamiltonian simplifies considerably: 

Hb — E -\- QnijiQn2j2 4{-^nijin2j2 ~ ^niiin2j2) 

~l~ \^QnijiQn2^nijin2 
+ 2-^n2ni{(ln2Qni Pn2Pni) 

and can even be written in a nice matrix form: 



(72) 



TT 7-1/ , / n \ ( 2 (^"iii"2j2 + ^nijin2i2) 1 / -Pn2j2 

H = E + -{Pn,n Pn,) [ 

Fn2ni I \ Pn2 



1 \ / 2 ^"^"iii"2i2 ■^n\jin2j2) ^'/^■Dnijin2 \ I Qn2j2 




, ^ I \ I 2 \ "iJi"2J2 - -ii.yjlii.2J2 J » i<.iji"2 1 I II.2J2 1 / Pro \ 

\ V '^-'-^n2ji2ni " n2n\ 

where E' is a constant that has to be chosen s.t. = Emf- This is a viable starting 

point for proving the equivalence between the gRPA formulation originating from the time- 
dependent variational principle and the one using the quasi-boson approximation. We will 
rewrite eq. ()73j) into a Hamiltonian that leads to the same eigenvalue equations as does 
the Hamiltonian of the small- fluctuation approach given in eq. ()A33|1 . From this we will 
conclude - since we have both times the same matrix to diagonalize - that the spectra of the 
Hamiltonians of the two approaches are identical, and we will see a correspondence in the 
eigenvectors. For this purpose we note the following: 



This excludes the cranking Hamiltonian Her = H — [ x oj'^ {x.)T°- (x.) - where H is the Yang-Mills Hamil- 
tonian and ^"(x) is a Lagrange multiplier, cf. e.g. [^I - here, since it is not of the form given in eq. ©. 
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1. We start with the term quadratic in p, P; there we have 

'2^^ + ^)nijin2j2 

1 1 1 1 1 1 

~ (2^^imi^jlii) (2^n2ni2^i2i2) (^■'^)™i«i;™2i2 C^^) 

and 

where we have used the abbreviation 

+ GiAr (76) 
This suggests that one should introduce new momentum operators 

Pmiii = Ij^riimJ-^jiii-^niji) (77) 
Pmi = ^f^nrmiP"!- (78) 

2. Now we consider the term quadratic in g, Q; there the observation 

^minj = ^UmkiUik2Unk3UjkiT7^ {T) (79) 



s 


6 


SGkik2 


SGk^^ki 


6 


6 


SGkik2 


5Gk,,ki 



Bminj ~ ^minj — ^UmkiUik2UnkzU jk^YT^ 77^ (^) (§0) 

(where T denotes the kinetic energy \'nf) is useful. Then one can rewrite 



5 5 

^UmkxUik2Unk-iUjkAY7^ 77^ {H). 

0Lrkik2 OLrkzki 



(81) 



Next, one can put to use the fact that the mean- field vacuum is a stationary point 
(s.p.) and thus 5{H)/5G = 0, the identity eq. ()E19|1 from app. lEl and the fact that 
the kinetic energy {■^irf) is independent of 0: 

-'^n2ni 2 "2".i 

= -^(T7^iT)]Uk,n,Uk2n2 
\0^kik2 / 

\0^kik2 J 



-{E) ) Uk,n.Uk2n2. (82) 
k2 
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This now suggests introducing new coordinates as well: 

q™, = V2Un^rn^qm- (84) 

3. One should note that the newly defined coordinates and momenta still fulfill the canon- 
ical commutation relations eq. (j71|) : thus, the new definitions amount to a canonical 
transformation that leaves the dynamics unchanged. 

We therefore end up with a Hamiltonian where the matrices in the quadratic forms are 
identical to those appearing in the Hamiltonian of the small-fiuctuation approach, eq. ()A33|) . 



H = E' + -(p P) ,1 (85) 




5G5il) 5G5G 

since H = {H) as defined in eqs. (IAin|l . ()A17j) . As a last point, we want to show that 
the eigenvalue equations that result from the Hamiltonian given in eq. (|85p are identical to 
eqs. ()A31|) . ()A32|) . This is in fact quite simple; since the excitation operator Q^^^ is a linear 
combination of B, B^, a, one can write it just as well as a linear combination of q, Q,p, P 
or q, Q, p, P. Thus, we may write 

Qbu = X!(<5LQmi + PLPmi) + ^{q'^flm + P^mPm) ■ (86) 
mi m 

We can now use the gRPA equations eq. (jlUjl to derive the eigenvalue equations: we insert 
eq. (jHSI) and eq. (|S^ into [HsyQ^'Bu] ~ ^^Q^Bu compare the coefficients of the various 
operators q, p, Q, P. The resulting equations are then easily expressed as 



S4>Srf> 5<j>5G 
&GS<t> SGSG 



'"^ '"-^ (Gl) 




(87) 



and 



54>54> 5(j)5G 





(88) 
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where Qs ^ Ps denote the symmetric part of Q'^, respectively, e.g. {Qs)ij = Qij + Q'^i = 
Qjjj}. These equations are now identical to eqs. (jXHTJ, ()A32j) . only the components of the 
vectors have acquired new names: 

6^^p'' ; 6G^ and Stt ^ q'' ; 5S ^ ^g,^ (89) 

Thus, we have proven that the two different approaches to the generalized RPA, namely the 
small-fluctuation approach from the time- dependent variational principle, and the operator 
approach, give in fact the same spectrum of possible excitations. Now the comments made 
at the end of app.^Jcarry over to the operator formulation of gRPA: If the stability matrix of 
the mean-field problem is positive, all 'eigenvalues' are larger than zero, and thus all fij^s 
are real; if we are not at a minimum of the energy with our choice of the mean-field vacuum, 
the stability matrix will also have negative eigenvalues, and thus we will obtain complex 
conjugate pairs ±2|fij/|. As long as we have real 'eigenvalues' fij, in the gRPA problem, the 
mean-field vacuum under consideration will be stable at least w.r.t. small fluctuations. One 
should note that this relation was of some importance in proving the bosonic commutation 
relations of the excitation operators. 



VI. THE MOMENT OF INERTIA 
A. General Form of Kinetic Energy 

From the discussion in sec. lV Al we saw that we can in general write the Hamiltonian in 
quasi-boson approximation as 

where {z^+} denotes the set of modes with positive f2y and {uq} denotes the set of zero modes. 
However, quite often an alternative expression to eq. (j9Up in terms of symmetry generators 
is useful. For this we have to recall from sec. llVl that if our reference state is a solution of the 
Rayleigh-Ritz equations, and if the symmetry generators under consideration are one-body 
operators, then the quasi-boson approximations P^ of the symmetry generators commute 
with the quasi-boson approximation of the Hamiltonian. This implies that they are zero 
mode solutions of the gRPA equations. If we now assume that there are no 'accidental 
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zero modes', i.e. the whole space of zero modes is spanned by the F^s, then we can write 
the PjyS, V G {z^o}, as hnear combinations of the F^s, and ^y(z^^y^^^Pl becomes a quadratic 
form in terms of the F^s. The advantage of this alternative expression is two-fold: First, in 
general the quasi-boson approximations of the symmetry operators are known. Second, the 
expression one obtains is (physically) more transparent, and can be more easily compared 
to other frameworks like e.g. the Thouless-Valatin method or the Kamlah expansion 



The general dependence of on the generators F^ will then be^^ 

1 



iiB,.m = ^r^M-T^'^, (91) 



where Ai is the moment- of-inertia tensor, and Einstein's summation convention is also used 
for the indices a, b which are employed to label the different symmetry generators F^. We 
cannot (as we have done in preceding discussions where masses were set to 1) 'normalize' 
Ai 'away', since the normalization of F^ is fixed by its very nature of a known operator. 
Furthermore, if the original symmetry generators (before the quasi-boson approximation) 
are generators of a non-Abelian symmetry, their normalization is fixed by the commutation 
relations. Thus, in this section, we will see how one can actually compute the moment-of- 

seen in sec. lVAI that in cases of appearances of zero modes 
one has to pass to a description of the oscillators in terms of canonical coordinates^^ and 
momenta (which we have done above). This actually allows to determine the moment-of- 
inertia tensor. We assume that we can construct a set of coordinates 9^, s.t. 

[Ql,r%]=zS^', (92) 

and which commute with all other normal modes. Since we already know the form of the 
part of Hb that is not supposed to commute with 9^, we also require^® 

[e%,HB] = tM'''Tl,. (93) 



HB,zm means 'that part of the Hamiltonian in quasi-boson approximation that only contains the zero 
mode operators'. 

Let us once again mention that the reduction of the non-Abehan to an Abehan symmetry (in the Yang- 
Mills case) in the quasi-boson approximation seems to be essential since otherwise one cannot construct 



canonically conjugate 'angle' operators fulfilling the canonical commutation relations [3£ 
We assume here that the moment-of-incrtia tensor is symmetric, which is at least true for the case of 
Yang-Mills theory, since there [r|, T^] = 0. 
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We start by making an ansatz of 0^ as a generalized one-body operator in quasi-boson 
approximation^^ : 

®_B ~ ^Q;miQmi + ^ P-mi^mi + ^q-m1m + ^p-mPm (94) 

whereas in general Tb reads 

^ B ^ Q;miQ mi "I" ^ P;mi^mi ~^ ^ g;m1m ~^ ^p-mPm- (95) 

In the latter case the coefficients are obviously known. We now insert eqs. (j94|) . (|95|) into 
eq. (HSl), and compare coefficients. This gives in general four equations, and allows to deter- 
mine 0Q.mi) 0p;mi5 0g;m5 ^p;m terms of A^^^ and the coefficients of F^. If we now insert 
the thus determined coefficients of into the equation that results from eq. we obtain 
an equation of the type 

{M-'T''^^'''' = i5''\ (96) 

where A/" is a matrix given entirely in terms of the matrices A, ...,F and the coefficients of 
r^. Thus, one has to invert the matrix A/" in order to obtain the correct kinetic term for 
the zero modes. One should note that the logic is just the other way around from what one 
would expect, namely that 0^ is defined via its commutator with Tb- However, for a proper 
definition of we also need the fact that it commutes with all the other normal modes. 
These general considerations find their application to the specific case of Yang-Mills theory 
in sees. IyTT!! IyTdI below. 



B. Energy Contributions of the Zero Modes 

We have discussed in some detail that zero modes cannot be treated via the ordinary 
creation/annihilation operator formalism; they have to be treated with the help of canonical 
coordinates and momenta. This changes their contribution to the gRPA vacuum energy. 
The gRPA energy in eq. ()43p was fixed s.t. 

{Hb) = Emf. (97) 
Since now Hb no longer contains oscillator modes only but rather has the form 

i^B = i?RPA + ^n(Ai-')"%+ ^"QlQ-^ (98) 

The usage of the p-q formulation as in section IV BI will be useful in this context. 
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where {z^+} denotes the set of modes with positive Qi,, eq. becomes in this context 

Erpa = Emf - ^{TUM-Y^b) (99) 

ue{v+} \ mi i 

which is again similar to the result obtained in nuclear physics In the special case of 
Yang-Mills theory to leading order in perturbation theory, we obtain an energy correction 
(see below) to the mean-field energy due to the zero modes^° 

Ai? = ^G2f(x,y)(r<^(x)r^(y)) (100) 

ihat is essentially identical in form to what one obtains in a second order Kamlah expansion 
^ and also to the Thouless- Valatin correction proposed in , but also with an important 
difference: the energy contribution which is due to the zero modes (and thus ultimately 
due to the deformation) is subtracted off after variation of the mean-field vacuum wave 
functional, and not before, as in the cases of P, 0, Is^; therefore, the determination of the 
parameters of the mean-field vacuum is not influenced by the correction. 



C. Special Considerations for Yang-Mills Theory 

We have repeatedly emphasized that the quasi-boson formulation works properly only if 
the reference state fulfills the Rayleigh-Ritz equations (the parameters are such that 'the 
energy functional is at its stationary point'). In Yang-Mills theory without a cranking term^^ 
this leads automatically to 

77 = and S = 0. (101) 
20 In this context, one should note that (r|r^) = {FT''). 

2^ As mentioned before the cranking Hamiltonian differs from the ordinary Yang-MiUs Hamiltonian by the 
term J cfi x uj"" {xjT'^ (x) , where a;"(x) is a Lagrange multipher. 
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This simplifies the expression for the Gauss law operator considerably, since, when we insert 
these results into the expression given in App. IC 21 we obtain^^ 



with 



pa _ pa p . _|_ pa „ 
^ B ~ ^ P,mi-'- mi ^ ^ p,mt'rn 

= (r?,(x))g(xi,X2)i^^(xi,X2) + (r^(x))^(xiK^(xO (102) 
X (f/-')x^/^"H^"')xfzf (103) 



where we don't integrate over x. Further remarks on notation can be found in footnote! 
One should note that all double indices are summed over except for x ! Inserting eq. ()94p 
into eq. ()93|) . with the Hamiltonian given in eq. (f?^ . we obtain (apart from Qp = <dp = 0) 



(105) 



Oa r)ll _ / A^ — IXafepb 
^Q.mi mini W ^ / 



where we have used the properties of the stationary point, i.e. the matrices A, F are all 
real, etc. [cf. sec. lV Bl and eqs. fl(^8j) - fl(n5j) ]. We invert the second line in a way that will be 
discussed below in some detail. Let us only mention at this point that (B^^)^^ is defined by 

jiniJ2»n-2 

2 ('^*i*2 '^"ii'"2 ~^ ^iim2^mii2^ ■ (106) 

We insert the result into the normalization eq. (jU^ and obtain 

{M-'r{n,m{F~')rmK. + 2T'p^^MB'')-'Un,T'p^n,) = 5"', (107) 



In this section we have to give up the super-index notation used until now in some places; instead of 
one super-index, the operators p will carry three indices (color, spatial, and position), i.e. Pm p\{p^) 
where b is the color, i the spatial, and x the position index, and correspondingly for Pmi- In this context, 
also the index a carried by has to be re-examined; in fact it is also a super-index, consisting of a 
color index a and a position index x: r^(x); the same applies to the super-indices that are carried 

by the moment-of-inertia tensor: Ai"-'^ A4"''(x, y). Having clarified this, in the more formal parts 
of this section, we will still stick to the super-index notation, since otherwise the formulas will become 
unreadable. 
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which bears some resemblance to the expression obtained in nuclear physics [16[, especially 
if one notes that {B^^)minj oc {A + B)minj and Fij oc {E + F)ij at the stationary point. 

D. Explicit Calculations 

Unfortunately, we cannot read off an explicit expression for {Ai^^Y^ from this in general. 
However, we can do two things: First, we can try to evaluate the terms 

^p,miP ^)mi^p^i and 2Tp,^^{(B^^) ^)minj^'p^nj (108) 

further. This will provide us with the moment of inertia - though for the energy correction 
due to zero modes we will need its inverse and this cannot be given in a general form. Second, 
we can go back to perturbation theory, where we in fact can give the inverse for the leading 
and next-to leading part in g"^ for the moment of inertia. 

It is a simple exercise to take ^ and {F~^)mi to calculate 

{M-'r{^,y) {{r;{y))l (zi)(F-i)^^^^4(zi,Z2)(rJ(z))^^^(z2)) = -(DD)^'=(z)(>i-^)='^(z,x). 

(109) 

The evaluation of this expression needs a bit more thought. The first step is to give a 
practical expression for {{B^^)^^)minj beyond its implicit definition in eq. p06|) . For this 
purpose, it is useful to rewrite G~^, cf. eq. (jT)), as 

A 

where are the eigenvalues of and P"^ are the projectors onto the corresponding 
eigenspaces. The projectors are complete in the sense that 

E P3 = (111) 

A 
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Then 

Using this expression for B^^ it is simple to verify that 

((-^ ) )iii2mim2 ~ / ^ ~ i \ C^nmi P«2m2 '^nm2 '^i2nii ) • (-'--'-4) 

In order to give a compact expression, using super- indices, for Tp^^- we introduce T^: with 
n,m being super-indices, where m stands for (xi,6, i) and n stands for (x2,c, j) , we define 
(no integration over x !) 

{f''T'S^,J^^2Si, = {f^'Un, (115) 

which inherits the property of anti-symmetry from {T"-)^^ = p"-'^ where p""^ are the SU{N) 
structure constants. Then we have 

{m^)Un2 = ^{{Uf^U^X,^, + (116) 

This expression can now be cast in the same form as (B^^)^^ by using the eigenvalue de- 
composition: 



namely one obtains 

(r?>(x))„,„, = ^-J2 ^^^(P^4"P'^)n2n.. (118) 
^ A,B ^^aXb 

Thus, we obtain overall 



2(Pp(x)),,,((i?ii)-i) 

«i«2m,im2 (rp(y))mim2 = -itE 



AaA_b(Aa + A_B 



-Tr(P^f^P^fJ) 



(119) 

Of course, the moment-of-inertia calculations ultimately have the goal to be connected to 
the static quark potential. It turns out to be interesting to compare the structure presented 
in eq. ()119|1 to the structures one obtains in a perturbative Coulomb gauge investigation into 
the static quark potential 
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3. A Short Excursion into Coulomb Gauge Perturbation Theory 



In Coulomb gauge, the charge density consists of two parts, namely the external part Pg^j 
and the gauge part p^auge — /"^'^A^(n*'')?, where II*'' denotes the transversal components of 
the momentum (the longitudinal component is eliminated using the Gauss law constraint). 
We then decompose the Coulomb gauge Hamiltonian into a free part, independent of the 
coupling, and a rest. In Coulomb gauge, the Hamiltonian contains arbitrary powers of g. 
The free part (i.e. that part of the Hamiltonian which is left if we set g — Gi) defines what in 
the following is to be called gluons. Since we are in Coulomb gauge these are given in terms 
of the transversal fields and the corresponding momenta. 

In order to compute the energy corrections to the energy of the perturbative ground state 
by the presence of external charges to 0{g^) in Rayleigh-Schrodingcr perturbation theory, 
we need on the one hand the 0{g^) contribution to the interaction Hamiltonian in order to 
compute 



and on the other hand the 0{g^) contribution to the interaction Hamiltonian in order to 

compute 



where the sum runs over all states except for the (perturbative) vacuum, and — -E^ is the 
energy of this state (the energy of the perturbative vacuum has been set to zero) . Since we 
are interested only in that part of the interaction energy that is proportional to Pexti^i) x 
Pea;t(x2), wc drop everything from Hint which is not quadratic in pext for the first term, and 
not linear in pg^t for the second. Since the term of the Hamiltonian which contains p at all 
is 



(which consequently has to be expanded in powers of g), for the second term we need 
Ptxti'^'^ + O^"-) p^g^^gg. Whereas the terms stemming from the first order contribution give 
the anti-screening part of the one-loop /9-function, the second order contribution gives the 
screening part of the one-loop /^-function (which is given by — ^ times the anti-screening 
part). We will go into a little bit more detail of this latter contribution. What goes really 



(0|i/i„t|0) 



(120) 




(121) 




(122) 
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into the computation of eq. ()12ip is 

= r'^miinfYiN). (123) 

Using the representation in terms of creation and annihilation operators it becomes clear 
that the matrix element {0\pgauge\^) is non-zero only if |A^) contains two gluons. If the 
gluons carry energy ujkiiOJk2 (and color index mi,m2) one obtains 

{N\pU9e\^) ^ n'"'''^^^y= (124) 

The comparison to eq. ()118|) is striking, especially if one takes into account that in the case of 
absence of condensates the eigenvalues of can be interpreted as single-particle energies, 
cf. app. |X| The second quantity needed is -Ej^ is obviously the sum of the energies 

of the two gluons, i.e. = ujki + ^^k2, and therefore the similarity of l/E^ to (B^^)^^ 
is obvious, cf. eq. Therefore, it appears that the "quantum part" of the moment of 

inertia contains only the (albeit generalized) screening component of the static interaction 
potential, and the (obviously dominant) anti-screening component has to be found elsewhere. 
This has to be subject to further investigations. 



4- Evaluation to Leading Order in Perturbation Theory 

After this excursion, we go back to perturbation theory to make an evaluation of our 
expressions for the two components needed for the moment of inertia. We note^^ that Tp is, 
in the perturbative scaling scheme given in app. IB 2\ of higher order in g than^^ Tp. Thus, 
we just compute the leading order piece of the moment of inertia: 

{M-T{^,y) (K(y))l;(zi)(F^^)^\''4(zi,z,)(rJ(z)):;^^(z2)) = 

-(DD)^^(z)(^-i)^'^(z,x) = (125) 

From this we conclude that the leading order piece of the moment of inertia in a perturbative 
expansion is just the Green's function of the covariant Laplacian in the background field A: 



{M-T (x,y) = -G'^(x,y) + Oig') (126) 



This can be read off the expressions eqs. (I1U3II . (|104|) together with the mean- field expectation value of 
r"(x) given in eq. HC16|) . 

For the purpose of counting powers of g, we take A to be of 0{g^^), s.t. D is completely of 0{g^). 
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with 

DSD5Gl^(x,y)=5-5.y. (127) 
Obviously, in the case of a vanishing background field, we obtain Coulomb's law: 

(>l-^)"^x,y) = f ^— . (128) 
47r |x — y| 

Let us now consider the perturbative evaluation of ea. pi9|) . If we solve the mean-field 
equation in a perturbative approximation (cf. e.g. [ll. IsilssI) for vanishing A, we obtain for 

(G-)^^(x,y) = 26^"! ^e^'^-(''-)|k|(5., - ^). (129) 

This shows immediately that we have a problem of building a gRPA treatment on top of the 
perturbative ground state, since has zero modes which is (as we have argued before) 
unacceptable for the gRPA formalism. Since this section has mainly illustrative purposes, 
we proceed by simply ignoring the zero modes in what follows. The non-zero eigenvalues of 
G^^ are given by A2|q| = 2|q|, the corresponding projector onto this eigenspace is given by 

(P^"^')S(x,y) = / ^-^^^nS., - ^)e^'^-(^-\ (130) 

where dfl^ denotes the solid angle in k-space. Inserting these expressions for A2|k| and 
(^p2|k|^aj, i^^Q (|TTT?|1 . we obtain 

2(r^(x)),,((5")-i),,^,„,(r^(y))„,„, = ^ln(A^^/A?^) J ^ye^^'^-y^S''^ 

+ finite terms (131) 

where finite terms are not determined explicitly. Am is an IR-cutoff needed in the calculation 
due to the primitive way used to evaluate the loop integral. The main point is here, however, 
that the "quantum part" of the moment of inertia is just the classical part, multiplied by 
a divergent constant. Therefore, abbreviating In (A^y /Af^) by a, we can now rewrite 
eq. (fTTlTjl as 

(-A,(A^-i)'^''(x,z))(l + a) = (5'^Xz, (132) 

where is the Laplacian w.r.t. the coordinates z. Thus, 

^"■^ 1 1 

{M-'r\^, z) = — - -. (133) 

^ ^ ^ ^ l + a47r|x-z| ^ ^ 
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If we introduce a renormalized coupling in the standard fashion jl| via 

1 1 



li vif 

+ |3o\n{^lJy/^^^) (134) 



the expression for Ai is made finite only for /5o = a, thereby corroborating our earlier claim 
that the moment of inertia can only account for the screening contribution of the static 
quark potential. 

It is interesting to note that the moment of inertia obtained here is the same as in a cranking 
type calculation j^. This also could have been anticipated, as it is well-known in nuclear 
physics that the moments of inertia determined from RPA and those from cranking are in 



fact identical in nuclear physics, cf. 



221. 



E. Interpretation of Moment of Inertia as Static Quark Potential 

We have to conclude that we cannot interpret the moment of inertia as the static quark 
potential directly. However, this is not necessarily fatal for a gRPA treatment of Yang-Mills 
theory. Due to the approximations made to the Gauss law operator (e.g. making it Abelian) 
it is not necessarily clear that^^ I7) = e*'^"®"|i?PA) is a correctly charged state, and thus 
(7|(A^~"'^)'*''r^r^|7) is also not necessarily the energy of a state with a certain prescribed 
charge. This demonstrates, however, that further developments are necessary in two areas: 
First, one needs a systematic expansion where the quasi-boson approximation is the leading 
order. Second, a clear interpretation of gRPA states in terms of real physical states appears 
necessary. 

The gRPA vacuum \RPA) is annihilated by aU gRPA annihilation operators Qg and by the gRPA 
approximated symmetry generators, Tb\RPA) = 0. Whereas in the case of non-zero modes the excited 
states are generated by the creation operators, in the case of zero modes the excited states I7) are 
given as plane waves and employ the canonically conjugate coordinate of the symmetry generator, I7) = 
e^''''^'' \RPA). These excited states are obviously eigenstates of the symmetry generators: r"|7) = 7°|7)- 
Their energy can be determined easily, too: Hb\j) = ^(A^~"^)"^7°7''|7). 
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F. Observations in Electrodynamics 



We want to close this section with a few remarks and observations on electrodynamics. 
In electrodynamics the (full) Gauss law operator is given by 

r(x) = vfn,(x). (135) 

We guess immediately the form of the corresponding "angle" coordinate as was discussed 
before in the context of gRPA; if we take 

e(x) = -Ga(x, y) Vf A,(y) with A-G'a(x, y) = (136) 

then on the one hand 

[e(x),r(y)]=^5.y (137) 

and on the other hand 

[e(x),iJ] = -G'A(x,y)r(y), (138) 
thereby identifying the moment of inertia as 

(A4-i)(x,y) = -GA(x,y). (139) 



It seems interesting that the angle coordinate eq. (I136D shown here appears to be closely 

nn 

related to the dressing function used in |4I|, cf. also 41], to construct gauge invariant 
electron fields, 

^jJ^{x.) = e'^^^("'y)^^^^(yV(x). (140) 

The Hamiltonian of electrodynamics furthermore allows for a decomposition very similar to 
the one made possible by the gRPA treatment for general bosonic theories, cf. eq. fj90|) . If 
we introduce the longitudinal and transversal projectors P^(x, y), P^(x, y) as 

Pf/x, y) = VrGA(x, y) VJ and Pj(x, y) = S,,5^y - Pg(x, y), (141) 

we can decompose the Hamiltonian into one depending solely on transversal degrees of 
freedom, nf(x) = P^(x, y)nj(y), Af(x) = P^(x, y)Aj(y) , and one depending solely on 
the Gauss law operator (symmetry generator): 

H=^jd^x (nnx)nf(x) + B2[A^]) +1 j rf3a;^d3^2G'A(x2,xi)r(xi)r(x2). (142) 
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The gRPA treatment, however, differs from the general treatment discussed above, though 
not in an unexpected manner. Since electrodynamics is eventually a free theory, the solution 
of the full Schrodinger equation gives the same result as the mean-field treatment, thus we 



can use as reference state the well-known ground state of electrodynamics 



ip[A] oc exp {-^ I d?xd?yAi{-K) 



(2^ 



p2 



)e 



ip-(x-y) 



A,(y) . (143) 



This state, however, is annihilated by the Gauss law operator and consequently we expect 
to vanish, and indeed this is verified by explicit computation. It should be noted in 
passing, however, that an application of eq. ()125|) for the moment of inertia gives the correct 
result (since the zeros in Fp, which come about since V acts upon the transversal U~^, are 
canceled by = G which is infinite in the longitudinal subspace, being proportional 
to the transversal projector). 



VII. SUMMARY AND CONCLUSION 



Let us shortly summarize what has been achieved in this paper, and which problems 
remain to be solved. We started by introducing shortly the canonical quantization and 
Schrodinger picture treatment of general bosonic theories with a standard kinetic term. We 
then considered the operator formulation of the (generalized) Random Phase Approximation 
that is the one common in nuclear physics, and which has recently also been investigated in 
the context of pion physics. We demonstrated that this approach can also be implemented 
for a general class of bosonic field theories, but that the class of excitation operators to be 
considered has to be larger than in nuclear physics; namely, one has to also allow for terms 
linear in creation/annihilation operators of the fundamental boson fields. It turned out to be 
possible (at least for a certain class of Hamiltonians with standard kinetic term) to prove the 
equivalence of the operator formulation to the formulation starting from the time-dependent 
variational principle. Then we demonstrated that, in the absence of zero modes, the gRPA 
Hamiltonian is just a collection of harmonic oscillators. The zero modes required special 
attention, but the problems could be solved along lines paralleling nuclear physics. We then 
considered the question of how conservation laws of the full theory translate into conser- 
vation laws of the theory in generalized Random Phase Approximation, and saw that, at 
least in the case of symmetries generated by one-body operators, existence of symmetries in 
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the full theory implies existence of symmetries in the gRPA-approximated theory, although 
these symmetries need not be the same; in Yang-Mills theory, the generalized RPA only 
carries an Abelian symmetry. This also opened up a possible way of improving the gRPA by 
requiring a certain fulfillment of the commutation relations by the approximated operators. 
Then we investigated the difference between the mean-field energy and the energy of the 
gRPA ground state with special emphasis on the energy contribution of the zero modes. For 
this purpose we had to calculate the moment-of-inertia tensor which (at least to lowest order 
in perturbation theory) turned out to be the Green's function of the covariant Laplacian in 
the background field at the stationary point. 

Furthermore, we compared the quantum corrections to the moment of inertia, and found 
that their structure is very similar to the screening contribution in a Coulomb gauge per- 
turbative calculation of the static quark potential. 

As a last point, we made some observations of electrodynamics where the concepts of mo- 
ment of inertia and zero mode operators can be applied to the system without any form of 
approximation. In the gRPA treatment of electrodynamics some of the general observations 
can be verified explicitly; the most interesting point seems to be that - even though the 
Gauss law operator vanishes in the gRPA treatment, due to the fact that the reference state 
is gauge invariant, - it is still possible to determine correctly the moment of inertia. 
To put the method into perspective, let us summarize the main positive and negative as- 
pects: 

On the positive side, we first have to mention that in the generalized RPA only energy 
differences w.r.t. the ground state are computed. This simplifies matters, since that part 
of renormalization that is usually done by normal-ordering is automatically taken care of. 
This brings us directly to the second point: energies of excited states can be computed. 
This is usually very difficult if one relies upon e.g. the Rayleigh-Ritz principle. The most 
important point, however, is the effective implementation of the Gauss law constraint. Even 
though the gauge symmetry is broken in the mean-field treatment, the unphysical excita- 
tions generated by the (gRPA approximated) Gauss law operator (more precisely given by 
plane waves given in the coordinate conjugate to the Gauss law operator) arc orthogonal to 
all the other physical excited states, which is almost as good as if they didn't even exist. 
But there are also a number of drawbacks of the method presented. The first has to do with 
the ground state energy: Whereas one does obtain a lowering of the ground state energy 
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w.r.t. the mean-field energy due to deformation, in a manner wliicli is even formally quite 
similar to the lowering obtained e.g. in the second order Kamlah expansion (at least to 
the order in perturbation theory considered) and also in other methods, the point of the 
calculation at which the energy correction due to the deformation is considered is funda- 
mentally different. In the Kamlah expansion and also in the Thouless-Valatin method, the 
corrections are subtracted before the variation is carried out, whereas in the generalized RPA 
they are subtracted only after the variation. Therefore the energy correction does not have 
any influence on the parameters of the mean-field. This brings us directly to the next prob- 
lematic point: of course the mean-field vacuum is not the gRPA ground state. In nuclear 
physics, one was able, however, to construct the gRPA ground state from the mean-field 
ground state; but, in the presence of zero modes, this state has a divergent norm. This is 
due to the fact that the gRPA is a 'small angle' approximation {2^, or in other words, the 
compact nature of the non-Abelian symmetry is lost. We have seen explicitly that in the 
case of Yang-Mills theory, the compact SU{N) symmetry is replaced by the non-compact 
?7(1)^^~^ symmetry. In nuclear physics, at least in the case of two-dimensional rotations, 
the fact that RPA is only a small-angle approximation was not so much of a problem, since 
the global dependence on the angle is generally known and therefore one can extract enough 
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information from the small-angle approximation to determine the whole wave function 
In Yang-Mills theory, we have no comparable knowledge that could be put to use practically 
and therefore, we cannot compute the ground-state wave functional. Lacking this knowl- 
edge, however, one needs other methods to evaluate matrix elements of operators, which 
up to now we haven't developed. Therefore, the only quantities we currently can calculate 
are the energies of the excited states. As a last though very important point, we have to 
mention that we have not dealt with the problems of renormalization. However, it should 
be possible to deal with them, since (at least in the case of 0^ theory) this problem has been 
faced already by Kerman et al Q| in the context of the generalized RPA derived from 
the time-dependent variational principle. 
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APPENDIX A: GENERALIZED RPA FROM THE TIME-DEPENDENT VARI- 
ATIONAL PRINCIPLE 



In this appendix we present the approach to the generahzed Random Phase Approxima- 
tion via the time-dependent variational principle as pioneered by Kerman et al 
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We show that it leads to equations of motion of coupled harmonic oscillators. The purpose 
of this appendix is to make this paper reasonably self-contained, so that the reader can 
follow our claim that (under the circumstances outlined in the main text) the operator ap- 
proach and the time-dependent approach to the generalized RPA yield identical equations 
of motion. 

We don't want to deal with a specific theory at the moment, we will only require the Hamil- 
tonian to be of the form 

H = \7:} + V[cl>il (Al) 

where (pi are a set of fields, and tTj are the canonical momenta conjugate to 0i, in the field 
(coordinate) representation under consideration 

I 64>i 

and i is a super-index, containing a position variable x, and all other indices required (like 
color, spatial etc). The Einstein summation convention is adopted, implying sums over all 
discrete and integrals over all continuous variables. V[(f)\ is a functional of the field operators, 
in the following referred to as 'potential'. 

The states that we consider as trial states for the time-dependent variational principle are the 
most general time-dependent Gaussian states (we only indicate the time-variable explicitly, 
i, j are super-indices and A/" is a normalization constant): 

^[<P, t] = ATexp (-{cP - 4>{t)u\G'\t) - ^S(t)).,(0 - 0(t)), + ^vr,(t)(0 - 0(t)),') . (A3) 
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The meaning of the parameters becomes clear by considering expectation values: 

imiM.im) = + G^Jii) (A5) 
= (A6) 



1 



(AT) 



We can now compute the action of the time-dependent variational principle 

S = j dt{i){t)\idt - H\i){t)) (A8) 

and obtain 

S= I dt< 

with 



J dtl TTi{t)^i{t)-ti{tG) + ^ti{G-^G) -nit)[(p,n,G,^\ (A9) 



n{t) = {ij{t)\H\m)- (AlO) 

We can now add a total time derivative^^ that does not change the equations of motion, and 
obtain for the action 

S = [ dt{T.,,{t)G,,{t) + irmUt) - nt)), (All) 



which shows that S is to be considered as the canonical momentum conjugate to G, and 
vf that of 0. The parameters of the wave functional are now determined via Hamilton's 
classical equations of motion: 

OTTi d(pi 

G.,{t) = ^ ; = (A13) 

However, in general it will be much too complicated to solve these equations, therefore we 
now consider a two-step procedure 

1. look for static solutions to the equations of motion 



Its form can e.g. be found in as (here only symbolically) ^(— | log (G ^) — EG). 
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2. consider small fluctuations around these static solutions. 



The static solutions are (obviously) determined via 

Ut) = ; ^i(t) = ; Gijit) = ; S,,(t) = 0. (AM) 

But these are nothing but the equations resulting from the Rayleigh-Ritz principle: 

— = 0- — = 0- ^ = 0- ^ = (A15) 
S-Ki ' S(j)i ' SEij ' 6Gij 

Thus, for a static solution of the time-dependent variational principle, the parameters are 
those which minimize (or at least extremize) the energy. This should not really come as a 
surprise, since for a static state the action reduces just to minus the energy times the 
respective time interval under consideration. 

For the next step, we decompose the general, time-dependent parameters into the static 
solution plus a time- dependent contribution that later on is considered to be small, e.g. for 
0: 

Mt) = ks + SMt), (A16) 

where (pi^g denotes the static solution^'', and S(f)i{t) the 'small' time-dependent part. We 
insert this decomposition into the equations of motion eqs. ()A12|) . ()A13|) . and obtain 

d-Ki 

= ^[0.,7r„a,S,] (A17) 

OTTjOTTi 

Now the meaning of etc. being small is clarified: in the equations of motion terms 
of higher than linear order are neglected (in the action, it would be terms of higher than 
quadratic order). The first contribution vf^, G^, S^] vanishes by virtue of the static 

equations of motion eqs. ()A14|) . ()A15|) . The same construction can be carried out for all four 
parameter types (0, tt, S, G), and the resulting equations of motion can be nicely summarized 



In other words, if we evaluate the first derivative 4?^ for tt = tTs, it is zero, and correspondingly for tj), G, S. 
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as follows: 



where ^ = 

wneie ^^^^^ SASB\ys' 
an auxiliary matrix Vl: 



f 6(P\ 

6G 
\5tJ 



( 





5^n 




5^H 




87^671 


SttSG 








s'^n 




5(f>54> 


54)S7t 


5(j)5G 












5Y;5(f> 


5T.5n 


SSSG 












5GS<i> 


SGStt 


5GSG 


5G5Y: 



J 



f 6(/)\ 
6G 



(A18) 



and can be made even more transparent if one introduces 



( 1 

-1 



1 

-1 oy 



(A19) 



Then one obtains 



/ 50 \ 

bG 
\bt) 



n 



\ 







s^n 




5<f>5<f> 




Scf>SG 


S<j)ST. 




S^H 


S^H 


S'^-H 






SttSG 


SttST, 






S'^H 


S'^H 


5GS(I> 


SGStt 


SG&G 


(5G5S 




S^H 


S^H 


s'^n 



J 



SG 



(A20) 



5S(50 <5S(57r <5S<5G ST,ST, 

At this point it becomes clear what determines the spectrum of small fluctuations around 
a static mean- field solution: it's the stability matrix of this static mean- field solution. Eqs. 
flA20|l are usually called the generalized RPA equations. Up to now, the only assumption 
that has been used was the assumption of being a Gaussian state. But we have also 
restricted the choice of Hamiltonians that we want to consider by eq. ()A1|) . This restriction 
will allow to carry the calculation a bit further. Ti depends only via (tt^) on both vf and S, 
V[(f)] depends only on (p and G. Thus for all Hamiltonians that we are considering, we have 
the kinetic energy written as 

. ,,1 



;7rf 1^) = + ^Tr(G-i) + 2Tr(EGE), 



(A21) 



where Tr denotes the trace over the super-indices. The static solutions vTj Sjj are deter- 
mined via 

(A22) 



and result in 



TT, 



0. 



(A23) 
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This information, together with the knowledge that V[(j)] does neither depend on vr nor on 
S, determines a number of second derivatives: 



5^H 



' 5Y;5(f> 



where we have used the abbreviation^^ 



^ 

^ 



■52?^ 

' §7Vi67Vj 

' <5S<5S 



(Gl), 



(A24) 



kl 



tr(SGS) 



— GkiSij + GkjSii + Gji6ki + Giidkj- 
These results can be used to simplify eq. (EJi 



(A25) 
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6n 
6G 
\6tj 
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1 
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SGS4) 







(Gl) 







(5G5G 



0/ 



6G 



(A26) 



with = 6ij. By taking the derivative of eq. ()A26|) with respect to time, and reinserting 



eq. ()A26|) . one obtains two sets of partially decoupled equations 
(I) . 

SG 




(11) 







Us i 





(Gl) 



SGS(I> SGSG 



&<j>5G 

{Gl) 




(A27) 



(A28) 



If we now make the ansatz of a harmonic time-dependence 

6G,j{t) ) \ 6G?^ 

sm 



(I) 
(11) 




(A29) 
(A30) 
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In this context, one should note that |^ — ^{SikSji + SuSjk) 
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we obtain the eigenvalue equations 

(I) UJ' 

(II) UJ^ 






5(j>54> S(j)&G 
SG54, 5G5G 




(A31) 



(A32) 



In analyzing properties of the generalized RPA equations eqs. flA31jl . ()A32|1 it is often simpler 
to study the Hamiltonian they are derived from than to study the equations themselves. In 
this context one can imagine eq. (fX2fi|l to originate from a Hamiltonian^^ H: 

This evidently is the Hamiltonian of a set of coupled oscillators. The important point is 
that the signs of the eigenvalues of the decoupled oscillators are determined by the reduced 
stability matrix'^^ containing only second derivatives w.r.t. to (f) and G (We assume here 
that (Gl) is positive definite and since {Gl) is only multiplied by objects symmetric in their 
two indices, e.g. 6T,ij = 6T,ji, this boils down to assuming that G is positive definite. This is 
a sensible assumption connected to the normalizabilty of ip[4>] and is discussed extensively 
in sec. nil a"]) . An interesting observation can be made immediately: usually will be 
only non-zero if there is a condensate (i.e. 7^ 0) in the system. Thus if we don't have a 
condensate, the equations for the one- and two-particle content^^ decouple, and, since in a 
system without condensate G will usually be translation invariant, we can see by considering 
the Fourier transformed quantities that l/g{p) just describes the energy spectrum of single- 
particle excitations with momentum p, where ^(p) is defined by the following procedure. 



They originate by the canonical equations of motion e p- At t = —SH/6{S(j>). 

The problem of coupled oscillators is well-known, cf. |43l l44l | . If the reduced stability matrix is a positive 
matrix (and thus our mean field vacuum is indeed a minimum), all the eigenvalues will be positive, and 
thus all oscillator frequencies will be real. 

We see in sec. lVBI that in a creation/annihilation operator formalism S(f>,5Tr are connected to the ampli- 
tudes of operators containing one creation/annihilation operator, whereas SG, S'S are connected to the 
amplitudes containing two creation/annihilation operators. 
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We introduce the Fourier transform G{p, q) by 

If we now require that G shall be translation invariant, i.e. (^(x, y) = Gijs. — y), we obtain 
for G{p, q) 

G(p,q) = (27r)35(p-q)^(p). (A35) 

This defines ^(p). 

APPENDIX B: HAMILTONIAN TREATMENT OF YANG-MILLS THEORY 
1. Yang-Mills Theory 

There is not much to be said about the Hamiltonian treatment of Yang-Mills theory in 
particular, since there is an excellent treatment in The main points to be kept in mind 
are: First, in order to work in the Hamiltonian formalism at all one chooses the Weyl gauge, 
i.e. Aq = 0. The price to be paid in this gauge is that the classical constraint equation 
D^^E^ = (E is the color-electric field, for the definition of D cf. eq. (IB4|1 below), i.e. 
the Gauss law equation, has to be implemented as a constraint on states: is a physical 
state if it is annihilated by the Gauss law operator (which is constructed from D"*E^ upon 
quantization). Even in the presence of this constraint the canonical pair of A, 11 can be 
quantized straightforwardly without the appearance of non-trivial metric tensors 46^ . 
To put our treatment of Yang-Mills theory in a nutshell: We consider it to be described as 
a canonical system, defined in terms of coordinates A°(x) and conjugate momenta n"(x) 
which satisfy ordinary commutation relations: 

[A1{^),U%y)]=t6^%,S^y. (Bl) 

The states in the physical subspace have to satisfy Gauss' law, i.e. 

r(x)|^) = 0, (B2) 

where we have introduced both the Gauss law operator 

r (x) = Df (x)n^(x) (B3) 
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and the covariant derivative in the adjoint representation 

Df (x) = ^'^^V. - {g)r'A'r{^). (B4) 

We have put the couphng constant in brackets since its appearance depends on whether 
we have chosen "perturbative" or "non-perturbative" scahng, cf. app. IB 21 The SU{N) 
structure constants are denoted as f"''"^. One should note that the Gauss law operator used 
here may differ from those in other publications by some proportionality factors, as it is 
not directly the generator of (small) gauge transformations. The reason for choosing this 
form of the Gauss law operator is that the only factors of g appearing during a change from 
perturbative to non-perturbative scaling appear inside the covariant derivative. 
The wave functional of the reference state taken for Yang-Mills theory is of Gaussian form 

.^[A] = Arexp|-(A,"(x)-A^(x)) l(G-i)^^(x, y) - zS^^(x, y) (A5(y) - A5(y)) | 

xexp{^7rf(x)(A^(x)-A^(x))}, (B5) 

falling into the class of reference states taken in the main text for generic bosonic theories. 
The remaining defintions, like the Hamiltonian, can be found in the next section, app. IB 21 



2. Factors of g 

In Yang-Mills theory one has basically two options concerning where one wants to put 
the coupling constant, either in front of the action (here called 'non-perturbative scaling'), 
or in front of the commutator term in the field strength (here called 'perturbative scaling'). 
In table lU we give a short list concerning which convention leads to which placing of factors 
of g in other quantities of interest. 





non-perturbative scaling 


perturbative scaling 


covariant derivative 




Df, = df,- igAf, 


field strength 


FjjLv — d^Ay — dyAfj^ — i\A^, Ay\ 
= i[D^,D,] 


Ff^u = df,Ay - diyAf, - zg[A^, A„] 


action 




S = -\F^^{x)F''^'^{x) 


electrical field 




E? = Fg, 
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non-perturbative scaling 


perturbative scaling 


magnetic field 


"R1 — — if:. ., pajfc 

= (V X A)^ - \P^%A'' X A*^), 


= (V X A)f - f /'^*'=(A^ X A^)i 


momenta Ilf = 




nf = i^-o = -E- 


Hamiltonian 


H = Sf-nf(x)nf(x) + :^Bf (x)Bf (x) 
= 2i^(E-(x)E-(x)+Bf(x)B^(x)) 


H = in"(x)nf (x) + |Bf (x)B^'(x) 
= i(E-(x)E-(x)+B-(x)B-(x)) 


wave functional of 
'free' theory 


^ A ~ e g 

' L J 


' L J 


generators of time-in- 
dependent gauge trafos 


[(-r«),(-r^)] = z<5,yr^"=(-r^) 


[|(-r^),|(-r^)] = |^'^xyr''^|(-r^) 


finite gauge trafos 
(gluonic part) 







TABLE I: Placement of the coupling constant g in the 'non- 
perturbative' and the 'perturbative' scaling scheme. We have 
used the shorter form for r"(x). 



APPENDIX C: EXPLICIT EXPRESSIONS 
1. Matrices A-F 

In this appendix, we want to give explicit expressions that are valid for the theories that 
have a Hamiltonian of the form eq. ©• As we have already mentioned, terms in the Hamil- 
tonian that contain more than four c/a operators do not contribute to any of the matrix 
elements due to the second gRPA approximation^^. Then the computation is straightfor- 
ward: we simply insert eq. (jltij) into eq. ()2(ij) . and compute the double commutators. The 

One can use Wick's theorem to compute the double commutators that appear in the definition of the 
matrices A, . . . ,F. Then one realizes that for those terms of the Hamiltonian that contain more than 
four c/a operators the evaluated double commutators still contain at least one c/a operator. Taking the 
vacuum expectation values - where according to the second gRPA approximation we use the mean-field 
vacuum state - of these expressions sets these terms then to zero. 
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computation is simplified by the observation that only those terms of H contribute to the 
matrices where the number of creation operators together with the number of creation op- 
erators in the definition of the matrix elements match the respective number of annihilation 
operators. This is the reason for the matrices A, C, D, E having contributions from one term 
of H only. B obtains contributions from two terms of H as is to be expected, whereas F 
obtains only one. Individually, the matrices read: 



Matrix A 



Matrix B 



= -24i/^r"'> (CI) 



with 



Matrix C 

r{jmi} 



Cjmi = -en'oi' (C4) 



Matrix D 



• Matrix E 



Matrix F 



Djmi — '^H^ ^ (C5) 



Er., = -2//o^f^ (C6) 



where ^'•'^ means: symmetrize in the indices i^j (i.e. add all permutations and divide by 
the number of permutations) and means 'that factor in the Hamiltonian that multiplies 
a creation and h annihilation operators'. We see that at the stationary point they simplify 
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considerably: 



B-minj = -^mmj "I" ^minj-i (C'S) 

22 ^"^ 

Bminj = 4:UmrmUii^UnniUjj^— — {V) , (Cll) 

OCPj^OLxrniii 

Djmi CjY^j, (C13) 



Emi = 0, (CM) 

1 

-( 

2 



-^mi — o^^i- (C15) 



2. Gauss Law Operator 

The Gauss law operator was defined in eq. ()B3|) . Here, we give its decomposition into 
creation/annihilation operators. After normal ordering, we obtain the result 



+ ^xiyi <^yizi l«zi +«zi , 



V SG"'^"'^'^^^^ \^-x.,i^^iyil J ^xizi '-^X2Z2 

( nlbi-ni ■\b2n2 I „feini 62^2 I 0„t^l"l/if'2»l2N 
^ I'^zi '^z, + '^Zl '^zo + ^'^Zl '^zo ) 



Here we have employed the notation 

^a^^aia2 _ j^aiaa2 (C17) 

where f"-^"-"-'^ denote the SU{N) structure constants. 
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APPENDIX D: COMMUTATION RELATIONS OF NORMAL MODES 



In this appendix we want to consider the question 'what are the conditions under which 
the eigenmodes can be treated as harmonic oscillators ?' We start with the observation 
that eq. (j4U|) implies by hermitian conjugation 

[Hb,QbA = -^vQb.. (D1) 

Thus, to every eigenfrequency, the negative eigenfrequency also belongs to the spectrum. 
Thus, without loss of generality^^ in the following we will assume that 

VL^ > 0, (D2) 

otherwise we just exchange the respective Qbu, Q'^bv "^^^ next point is that all the commu- 
tators 

[Q\.,Q\,] [Qbu,Q\,] [Qb.^Qb,] (D3) 

are pure numbers. This is due to the fact that all operators considered in this context are 
by construction linear in B, B"'', a, a) . We denote these numbers as 

[Q\,^Q\J\=M,.; [QB,,Q\^=N^,■ [Qb,.Qbu\=0,,. (D4) 

Up to now these numbers are arbitrary; one can, however, put the equations of motion to 
some good use. Consider 

^AQb,,Q^b.] '"'^ [QbAHb.Q^b.W (D5) 

JacoU ^^t^^^ ^^^^j _ j^t^^^ ^^^^^ ^^^^ ^j^g^ 

= [£BjV]+[gL,[£B^QBg] (D7) 

= ^,[Qb,^Q^bu]- (D8) 

In other words: 

{Vty — VL^)N^y = (no sum over /i, z/). (D9) 



■^•^ We have seen in sec. lV BI that the assumption that the eigenenergies of the modes are real aheady implies 
that we are deahng with a stable mean-field solution. 
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If all eigenvalues are distinct and non-zero then eq. ()D9|) implies that Nfj^^, is diagonaL Since 
the gRPA equations are homogeneous equations, we may now normahze the amphtudes in 
such a way that^^'^^ 

N^u = V- (DIO) 
The same procedure can be carried out for M^^^, O^^. We derive in analogy to eq. ()D9|) : 

{Q, + n^)M^, = 0, {n, + n^)o^, = o, (dii) 

where again no sum over double indices is performed. From this we conclude, again if there 
are no zero modes, that 

M^, = 0, O^, = 0. (D12) 

Thus, under the aforementioned conditions, the normal modes satisfy ordinary c/a commu- 
tation relations. To arrive at this, we have indeed used both gRPA approximations, since 
the basic ingredient was that is a one-particle operator (first gRPA approximation) and 
that all the commutators of a, a\ B, B''" are pure numbers (quasi-boson approximation). 



In it is stated that one can show that N^^i, is positive if one is considering a positive definite Hessian 
at the stationary point of the mean-field problem; a possible proof of this statement works as follows: 
only if the Hessian is a positive matrix it is guaranteed that all eigenvalues fl'^ are positive. This has 
been assumed so far (e.g. how we concluded that for every positive frequency there is also a negative 
one etc.); thus we know that the system is stable. Knowing this, we can argue as follows: we know 
that [H,Q''gJ = il.^Q^g^,[H,QBi,] = -^uQbv and [QBu.Q^By'] = J^u^^v' ■ From this we conclude that 
the Hamiltonian has to look like H — '^^{^v/Mv)Q^gj^QBv Now we can study two different scenarios, 
namely can be positive (we call the set of u for which this is true v+) or negative (correspondingly 
V-). In order to have the usual creation/annihilation commutation relations we have - for v ^ - to 
interchange 'creation' and 'annihilation' operators; for clarity, we introduce new letters for them, i.e. for 
V G v-,Qbv Pbv^Q\v ^ Pbv Now we normalize Qb^ and Pes s.t. for v £ i^+,A/l/ — 1 and for 
ly e iy-,J\f,y = -1. The Hamiltonian then reads H = J2„e„_^_ ^i^Qb uQb + J2iyei^-(~^'')PB i^^B ,y + const. 
But here we see that we are dealing with a rather unstable system: the more modes are generated by P^, 
the lower the energy becomes. This cannot be true, however, since we know that we started from a stable 
system (with all frequencies real) ! Thus has to be empty. 

This relation has also as its consequence that the states generated by Q^^, Q^bu from the gRPA vacuum 
(defined to be annihilated by all gRPA annihilation operators, Qb ij\RPA) = OViy) are indeed orthogonal: 
(RPA|QbmQL|RPA) = (RPA|[Qb^,QL]|RPA) = {\[Qb^.,QIJ\) = S^, where we have used the second 
gRPA approximation in the next to last step . 
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APPENDIX E: GENERAL POTENTIAL 



In this appendix we will demonstrate that a potential that is an arbitrary polynomial in 
the field operators, can - once one decomposes the field operator into creation and annihi- 
lation operators - be written in normal ordered form s.t. all the coefficients appearing in 
front of the normal ordered products of creation/annihilation operators can be written as 
functional derivatives w.r.t. G and (f) of the vacuum expectation value of the potential (the 
creation/annihilation operators are defined with respect to that vacuum). 
For simplicity we take the potential of the form 

=7W^,..,^„0^, (El) 

where Xi are super-indices. It is clear that, if the claim holds for this potential, it will hold 
for an arbitrary polynomial since it will be a sum of terms of type fjEl|) . 
Since the field operators commute, it is sufficient to consider an A^xi...a„ that is symmetric 
in all indices. We have seen in section IIII Al that one can write 

(i>.=4>X + U^y{hl + hy), (E2) 

where x, y are super-indices, U'^ = G, and = S^y For our purposes, it will be more 

useful to define rescaled operators 

al = U^yhl ; = U^yby with [a^, 0^] = G^y. (E3) 

Sometimes we find it also useful to write 

4>x = 4>x + ^x with Lp^ = a}^^ a^. (E4) 

After all this notational introduction, let's come to the proof. We note that, similar to 0, 
both and ip commute; thus we can write equally^^ instead of eq. ()E1|) : 

V{(\)\ = M^,,„^„ ^ 0X1 ■ (Pxrufx^ + l ■■■fx^- (E5) 

m=o\m J 

It is understood, obviously, that the index of x increases from left to right; if it should ever decrease, as 
in the case 771 = 0, the (f)S are to be considered absent. 
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Now we normal order the terms with a fixed m: 



1 + at ) ■ ■ ■ (a^ + at ) 



+ ax™+iax„+2---ax„. (E6) 

We see that every row contains a fixed number of creation/annihilation operators and that 
it contains all possible permutations of types (creation or annihilation) among the possible 
indices. Especially for every term 

■ ■ ■ «xp ■ ■ • flx, • • ■ (E7) 



(E8) 



there also exists a term 



with all undenoted operators identical. We now use Wick's theorem, cf. e.g. jlG], to put 
every line into normal order. It is practical to deal with the whole line for the following 
reason: if we use Wick's theorem naively we obtain the normal ordered expression plus the 
normal ordered expression of two less operators times their contraction etc. . However, a lot 
of these contractions are zero, since (we denote the contraction of two operators by C) 

C{a\a,) = 0. (E9) 

However, if we deal with the complete line at once, we can use that, since for every ar- 
rangement eq. ()E7|) there also exists a partner eq. ()E8|) . we will always obtain contractions 



C{alaj + ttia]) = dj. (ElO) 

It is clear that for k contractions present we need 2^ partners to obtain a non-vanishing 
contribution. The important point is that they exist if we deal with the whole line at 
once. Remember that the contractions are multiplied by normal ordered terms, and that 
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[a, a] = [a^^ , = 0. Since in addition A4xi...x„ is symmetric in all its indices, all terms with 
a fixed number of contractions, creation, and annihilation operators, will give an identical 
contribution. Thus the question appears: assume we start out from the line where every 
term contains p creation operators, q annihilation operators and consider now terms with k 
contractions; how many terms will we obtain ? The answer is simple, 

number of terms = — r x — — x I — — — - , (Ell) 



2^ p\q\ \k\ {p - ky. {q - ky. 

and comes about as follows:^^ 



1. the line containing p c's and q a's contains ^^{^ terms; 

2. out of the p c's and q a's we take k each, and put them together to form contractions: 
there are obviously ^^^'^^i (g^fc)i ways to do this. 

3. However, it does not matter in which order we perform the contractions since the 
contractions commute; thus in the step before we have overcounted by a factor of k\. 

4. We now have to take into account what was said above: On the one hand, a lot of 
contractions are zero. On the other hand we can form pairs - this gives as argued 
above an additional factor of 2~''. 

With this formula at hand, we can at first answer the following important question: assume 
that we have started from an expression with n field operators; upon normal ordering we 
obtain expressions with n,n — 2, . . . ,n — 2k c/a operators each; the question now is: is the 
relative number of P c's and Q a's (with P + Q fixed) always the same, no matter from 
which n one starts and how many contractions one needs^^ (provided n — {P -\- Q) is even) 
? This question can be answered in the affirmative in the following way: We start out 
with an expression that contains p c's and q a's; after k contractions we will end up with 
P = p — c's and Q = q — k a's. Since in the beginning p + q = n was fixed, and we end 
up with P + Q = n' fixed we need for every term the same number of contractions, namely 
2k = n — n' . With this we can rewrite eq. (jEllll as 

n\ 1 n\ 1 



2''k\ P\{n' - Py 2^k\P\Q\ 



(E12) 



•^^ In the following, we abbreviate 'creation operators' as c's and 'annihilation operators' as a's. 

As an example: no matter where one starts - if one can get to the expression with, in total, two c's and 
a's, will they always come as aa + a^a^ + 2a^a ? 
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Thus, we have decomposed the number of terms into a factor that depends on n, n' which is 
a constant for P + Q = n' fixed and P varying, and a factor that depends on the number of 
creation and annihilation operators. If we go back to eq. ()E5|) we see that we have different 
possibihties to end up with P c's and Q a's (P + Q = n'): either start from n ips, and 
perform k contractions, or start from two 0s, {n — 2)(fs and perform k — 1 contractions etc. 
Thus the contribution to Mxi...x„4>xi ■ ■ '4'x„ containing n' c's and a's can be written as 



Mx,...x^ 



1 

— TT rrrraL ■ ■ ■ ol a 

P\r,'- 



.P=0 



P\{n' -P)\ 



t ...„t 

XI "'Xp^'Xp+i 



^ ^ \2'^k\{n - n' - 2k)\ '"^ 

where the latter sum runs to — n') ii n — n' is even, and to j{n — n' — 1) if it is odd - 
we will deal with these details below. This expression allows us to write down the vacuum 
expectation value (VEV) of ^[0] since it corresponds to the case n' = 0. Distinguish 

• n = 2N: the VEV reads 

N , 



fc=0 

n = 2N + 1: the VEV reads 

N 



M\{n-2k)\ 



fc=0 

G 



2^k\{n - 2k)\ 



For the following treatment, we can treat both cases with the same formula if we realize 

that ■M.xi...x2N+i'i'x2N+i same properties as M.xi...x2n ^^"^ ^^^^ ^1- ^-b]15|) 

always at least one survives. Thus we only have to treat the case with n = 2N. We now 



•^^ Note that the factors of the contractions and of the binomial decomposition of eq. (jE5p can be put together 
in a practical manner: — 
(f> opposite to eq. I|E5|I . 



in a practical manner: ^"atfcf'" L+i-( J+2fc+i)) = 2fcfc!(n "ri'-2fc)! ■ ^1^° ^^at we have arranged here 
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consider 



5G 



iv[<p\) 



y\y2 
5 



5G 



2/1 i/2 



N , 
fc=0 ^ 



2^k\{n ~ v! - 2k)\ 



GxiX2 ' ' ' ^X,2k-1)^ 



n'=0 



{2k-l)^2kr'X2k+l r'X2N 



{*) 



M 



yiy2Xz...x2N 



N 



^ \ 2m{n - n' - 2k)\ 



n'=0 



k Gx^X4 ■ ■ ■ G X(^2k-l)X2k^X2k+l ' ' ' 4'X2N 



27V' ■^yiy2x„i+-^-..x2N 



X 



N~N' 

E 

fc=0 



n! 



2''k\{n - n' - 2k)\ 



G 



■■■G 



]^^„' + l^„'+2 ^ X(^2k-l)X2ky^X2k+l 'rX2N 



n'=2=2N' 



(E16) 



In (*) we have used the symmetry of Ai. Obviously, as has been indicated by the suggestive 
notation, one is not restricted to one functional derivative but one can also perform N' of 
them, and then the restriction A^' = 1 in the last line is rendered unnecessary. We see clearly 
that, apart from the factor 2^^' the outcome of A^' derivatives of the vacuum expectation 
value w.r.t. G is identical to the prefactor of the addend containing P + Q = n' = 2N' c's 
and a's in eq. ()E13|) . Thus we can rewrite eq. ()E13|) as 

5 5 



y I 



4i ■ ■ ■ ^yp^yp+i " '^yn 



,p=o 



^^yiy-i ^^yn'-iyn' 



{Vm. (E17) 



The treatment we have presented up to here is valid for P + Q even. If P + Q is odd, we have 
to perform derivatives w.r.t. (f) and thus we have to treat the n even/odd cases individually. 
Let's start with n even: 
6 



■{vm 



k=0 ^ 



nl 



2>'k\{n - n' - 2k)\ 



n'=0 



GxiX2 ' ' ' ^X(^2k-l)X2k^X2k + l ' ' ' 0x„_i0x„ 



(*) 



Mx,...X„ 



-12/1 

fc=0 
N-l 



2^k\{n -n' - 2k)\ 



^X\X2 ^ 



n'=0 



X(2k-V)X2k^ X2k-y 



Mx^ 
M 



yiX2...x„ 



12/1 

fc=0 

N-l 



nl 



k=0 



2''k\{n -n' - 2k)\ 



GxiX2 ■ ■ - G 



x(2k-i)X2k^X2k+\ 'rx„-i 



n'=l 



2f'k\{n - n' - 2k)\ 



r' c 



+ {2fc-l)'^n' + 2fc "'^n'+2fc + l 



n'=l/ 

(E18) 
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The main point happened in hne (*) where the fact that we started from even n played a 
role. If n is even, k — N means that this addend doesn't contain a single factor of 0, thus 
its derivative vanishes. This is different in case of n odd, there the upper boundary is not 
affected by the first differentiation. It is different if one performs two derivatives w.r.t. 
since then the upper limit of the sum changes once altogether independent of whether one 
starts from n even or odd. Thus one obtains a relation between derivatives w.r.t. G and to 

^ -WD = 77^ WD, (E19) 



2 6(l)xS(f)y 6Gxy 
which will be very useful in proving the equivalence between the generalized RPA using 

the operator approach and generalized RPA as derived from the time-dependent variational 

principle. 

To put this appendix in a nutshell, we have proved that a potential that is an arbitrary 
polynomial in the field operators can be decomposed into creation and annihilation operators, 
s.t. upon normal ordering one obtains a sum of subsums where each subsum contains a fixed 
number of c's and a's. The subsum consisting of the terms containing n' c's and a's can be 
written as a standard polynomial in c's and a's 



n ^ 



P\(n' - P) 
.P=0 ^ ' 



-at, • • • at„aj 



(E20) 



multiplied by 

• if n' is even, n'/2 derivatives w.r.t. G times a factor 2"'/^ 

• if n' is odd, one derivative w.r.t. and {n' — l)/2 derivatives w.r.t. G times a factor 

2(n'-l)/2_ 

We have also shown that each derivative w.r.t. G may be traded for two derivatives w.r.t. 
0. 
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